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ABSTRACT 


The phenomenon of cavity flow oscillation is investigated to determine the 
conditions for onset of periodic oscillations and to understand the relationship 
between the state of the shear layer and the cavity drag. Experiments have 
been performed in a water tunnel using a 4" axisymmetric cavity model instru- 
mented with a strip heater on the nose cone and pressure taps in and around 
the cavity. A complete set of measurements of oscillation phase, amplitude 
amplification along the flow direction, distribution of shear stress and other 
momentum flux is obtained by means of a laser Doppler velocimeter. Drag 
measurements were made by integrating the mean pressure over the solid sur- 
faces of the cavity. Results indicated exponential cavity drag dependence on the 
length of the cavity. A jump in the cavity drag coefficient is observed as the cav- 
ity flow shows a bluff body wake type behavior. An independent estimate of the 
drag, which is obtained by integration of shear and mean momentum transfer 
terms over the peripheral area of the cavity, confirms the exponential depen- 
dence of drag on the length of the cavity. Results, also reveal that the drag of 
the cavity in the non-oscillating mode is less than the case if the cavity were 
replaced by a solid surface. Natural and forced oscillations of the cavity shear 
layer spanning the gap are studied. The forced oscillations Eire introduced by a 
sinusoidally heated thin-film strip which excites the Tollmein-Schlichting waves 
in the boundary layer upstream of the gap. For a sufficiently large gap, self- 
sustained periodic oscillations are observed,while for smaller gaps, which do not 
oscillate naturally, periodic oscillations can be obtained by external forcing 
through the strip heater. In the latter case resonance is observed whenever the 

forcing frequency satisfies phase criterion = N, and amplitude exceeds 


/ 
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certaln threshold levels, but the phenomenon, is non-self -supporting. The drag 
of the cavity can be increased by one order of magnitude in the non-oscillating 
case through external forcing. For naturally occurring oscillations, it is possible 
for two waves to co-exlst in. the shear layer (natural and forced). Also, it is pos- 
sible to completely eliminate mode switching by applying external forcing. For 
the first time a test is performed to cancel or dampen the amplitude of Kelvin- 
Helmholtz wave in the cavity shear layer. This is done through introducing an 
external perturbation with the same frequency of the natural component but 
having a different phase. Seduction by a factor of 2 is obtained in the amplitude 
of the oscillation. 
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Chapter 1 

Introduction 


1.1. Background 

Cavity flow phenomena occur in a variety of structures; slots between mov- 
able parts of control surfaces of aircraft, ships, tandem arrangements of bluff 
bodies such as adjacent tall buildings or tractor-trailer combinations, gas 
dynamic lasers, hydraulic gates, and control valves. The presence of cavities can 
give rise to drastic changes in the steady and unsteady loading character of 
these structures. The flow past cavity-type geometries (cut-outs) may introduce 
large amplitude velocity and pressure fluctuations in and around the cavity 
which are often damaging to the structures. Flow oscillation over gaps between 
insulation tiles used on the space shuttle outer surface, caused by evaporation 
of the filler material at reentry, is suspected of being responsible for the loss of 
a portion of the tiles. Emission of acoustical noise is another unwanted aspect 
of cavity flow-type oscillation. On the other hand, flow over cavities can be 
effectively utilized to alter heat transfer (White 1971) and drag characteristics of 
aerodynamic configurations (Koenig 1978). 

Charwat, et al. (1961) categorized cavity flows into open or closed types : 
"open" when the separated shear layer spans the gap, while "closed" refers to 
flows where the separated shear layer reattaches on the bottom of the cavity 
and separates again ahead of the downstream corner. For supersonic speeds 
and for subsonic speeds with turbulent boundary layers, closed cavities are 
found to have a width-to-depth b/h ratio of more than 11 (Charwart, et al. 
1961). Sarohia (1975) shows a demarcation line between an open and closed 
cavity to be around b/h & 7-8 for low subsonic speeds with laminar upstream 
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boundary layers. Furthermore, open cavity flows are divided into shallow and 
deep cavities, depending on the type of acoustic resonance present in the flow. 
When an acoustic resonance is present in the flow, the standing acoustic wave is 
oriented longitudinally for shallow cavities (b/h > 1). The standing wave is of 
transverse type for deep cavities (b/h < 1) (Charwat et al., 1961, Plumblee et al., 
1962, East 1966, Tam and Block 1978). Such a definition for the shallow and 
deep cavities loses its meaning in the absence of the acoustic resonance in the 
cavity. Tani et al., (1961) define cavity flows as deep and shallow cavities based 
on their respective low and high pressure drag values. The demarcation line 
based on Tani's definition is about b/h » 1.4. 

Sketch 1 illustrates the physically significant parameters of the cavity flow. 



where U e is the upstream edge velocity outside the boundary layer, <5* 0 is the 
displacement thickness and 6 0 is momentum thickness. For fixed flow parame- 
ters (U 9w 9 0 , 5* 0 ), cavity flow shows a strong oscillation of the shear layer when 
the width b exceeds a critical minimum width, bmin> (Karamcheti 1956, 
Sarohia 1975, Rockwell and Naudascher 1978). The existence of this minimum 
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width for the onset of oscillation is linked to the nature of the instability growth 
in the shear layer. Previous investigators noticed that the main driving force of 
the acoustic resonance condition in the cavity was the amplification of instabili- 
ties by the respective cavity shear layer. In the shear layers, the final growth 
and saturation of the instability waves result in roll-up of the shear layer and 
appearance of large vortical structures (Brown and Roshko 197,4). The velocity 
fluctuation associated with the cavity flow oscillations can be attributed to the 
presence of these large vortical structures (Sarohia 1975, Rockwell and Knisely 
1979). In his study of the frequency selection process, Sarohia (1975) concluded 
that the frequency at which the cavity chooses to oscillate is such that: 

(1) certain phase criteria necessary for the self-sustainment of the flow are 

satisfied: and 

(2) maximum integrated amplification through the cavity shear layer occurs. 


The cavity flow oscillation was described as a self-sustained type. Rockwell and 
Naudascher (1978) suggest the following chain of events to be responsible for 
the self-sustaining nature of the oscillation: impingement of organized vortical 
structures upon the downstream corner; resultant upstream influence (inter- 
preted as Biot-Savart induction): reception and conversion of this upstream 
influence into velocity perturbation by the sensitive region of the shear layer 
close to the upstream corner; amplification of these fluctuations by the shear 
layer; growth into vortex structures. Consideration of the acoustic speed 
becomes important in the nature of the upstream propagating waves for high 
Mach number flows or when the impingement length, b, is of the same order of 
magnitude as the acoustic wave length (ft » \a) (Heller and Bliss 1975, Rossiter 
1964, and Rockwell and Naudascher 1978). 
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The self-sustained nature of the cavity flow oscillation is similar to that of the 
jet-edge oscillation (edge tone phenomenon). Self-sustained oscillation of cavity 
flow has the following properties: 

(1) for a fixed cavity depth and fixed flow parameters ( C/ e ,$) $o) there exists a 
minimum width (6mm) below which no oscillation occurs; 

(2) for fixed cavity depth (h) and fixed boundary layer thickness there exists 
a minimum Mach number or free stream velocity ( U e ) below which no oscillation 
occurs ; 


(3) for 6>6 min or{7 e > f/ e and fixed boundary layer thickness the fre- 

U g 

quency of oscillation decreases as the parameter — — decreases. 

o 


At certain values of the frequency undergoes a jump to a higher mode 

of oscillation with a corresponding higher frequency. Hysteresis is observed in 

U„ 

the transition region. In the new mode, frequency also decreases as -g~ 

decreases. Sarohia (1975) has found that for fixed values of U e , b, and h, there 
exists a maximum boundary layer thickness or (5 0 *max] above which no 

oscillation occurs. 


The ladder-type behavior of frequency versus 


Hi 

b 


(which has the dimension 


of frequency) is remarkably similar for self-sustained flows and externally 
forced free shear layers (Ho and Huang 1982). However, there is a significant 
difference between the nature of perturbations introduced at the sensitive 
region of the shear layer by the downstream corner and that of an external 
forcing source. External forcing subjects the shear layer to a constant ampli- 
tude perturbation at one frequency. But flows excited by the impingement of 


-5- 


the shear layer at the corner are subjected to variation in strength of events at 
the corner (Knisely 1980). Often this variation gives rise to the amplitude modu- 
lation of the fluctuations, which shows as a low frequency component in the 
spectrum of fluctuations. Knisely (1980) shows that transverse positioning of 
the successive vortices as they approach the downstream corner yields vary- 
ing incident structure severing degrees. It is not clear what sort of mechan- 
ism is responsible for the presence of these low frequencies in the cavity shear 
layer. Also, the precise role of the downstream corner as a source of external 
forcing for the shear layer and generation of different frequency components of 
oscillation is not fully understood. 

As mentioned before, the presence of cavities (gaps) can be used effectively to 
change the drag characteristics of bluff bodies. A novel experiment by Roshko 
and Koenig (1978) demonstrated that by properly choosing dimensions, a gap 
generated by placing a disk in front of a circular cylinder can drastically reduce 
the drag coefficient of the system. A similar observation for different 
configurations of bluff bodies in tandem was made by Mair (1965) and Zdravko- 
vich and Pridden(1977) . It is unclear into which category these oscilla- 
tions fall. Scattered information about the frequency of the shear layer oscilla- 
tion shows behavior characteristic of self-sustained flow for the low drag cases 
of flow over gaps generated by bluff bodies in tandem. 

1.2. Statement of the Problem 

A large portion of cavity flow investigations are devoted to the study of insta- 
bility aspects of this problem. Quite independently, a second group of studies 
can be found which deal only with the pressure field and drag characteristics of 
this flow (Roshko 1955, Tani 1961, Fox 1964, Koenig 1978). To our knowledge no 
investigator has tried to fill the gap between the first and the second group of 
experiments. It was this gap that motivated the present investigation into the 
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role of instabilities in drag and pressure field formation of cavities. 

Another important aspect of cavity flow study is to develop a technique to 
control the frequency and amplitude of oscillations. Previous work involving 
control of the cavity flow oscillation has dealt with changing the geometry of the 
downstream or upstream corner of the cavity by rounding the corners (Ethem- 
babaoglu 1973) or by placement of external objects in the cavity such as flaps, 
deflectors or perforated plates (Heller and Bliss 1975 and Ball 1957). Usually 
these devices introduce their own complexities into the problem. There was a 
need for a clean technique that can control the oscillations by purely introduc- 
ing external perturbation to the shear layer. 

The following were objectives set for the present investigation. 

(1) To understand the relation between the state of the shear layer and the 
drag of the cavity. 

(2) To determine the influence of initial conditions on the onset of self- 
sustained oscillation and to obtain a better understanding of the feed-back 
mechanism. 

(3) To develop and use a technique which is capable of controlling (enhancing 
or damping) established oscillations. 

1.3. Approach 

The connection between the state of the shear layer and the drag of the cav- 
ity system can be made by studying the quantities that are responsible for the 
transfer of the momentum in and out of the cavity (tTir andiTu). These quanti- 
ties carry the necessary information about the nature of the instabilities in the 
shear layer and play an inherent role in the formation of the pressure distribu- 
tion inside the cavity system. Thus, accurate measurements of these quantities 



for different width-to-height ratios (b/h) in conjunction with the pressure dis- 
tribution measurements would provide the necessary information for fu lfillin g 
the first objective. External forcing by strip heating was used to study the feed- 
back mechanism and control of cavity flow oscillation .thus, fulfilling the second 
and third objectives. 
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Chapter 2 

Experimental Details 


2.1. Water Channel 

The experiments were performed in a newly constructed free surface water 
tunnel, which is located at the Jet Propulsion Laboratory, California Institute of 
Technology. This tunnel is specially designed to attain a low turbulence level for 
the purpose of carrying out sensitive instability-type experiments. The essential 
features of the free surface tunnel are shown in figure 1. Water is recirculated 
in a closed loop by a 50 horsepower axial flow pump, located 5m below the free 
surface, in the lower leg of the tunnel. The desired flow speeds Eire achieved by 
throttling flow through an on-line and a by-pass valve. The upper portion of the 
tunnel consists of a diffuser, a straight wall section comprised of flow manipula- 
tors and a test section. The high speed flow from the pump is then diffused by 
the 22° half-angle diffuser. In order to prevent severe flow separation, the effect 
of an unfavorable pressure gradient is reduced by matched pressure drops of 
three perforated plates installed in the diffuser section. Before entering the test 
section, the flow passes through different flow manipulators and a final 4:1 fifth- 
order polynomial contraction. Each manipulator and contraction is designed to 
break up and suppress different scales of the turbulence. 

The test section is 45.7 cm (18 inches) wide and 63.5 cm (25 inches) deep. 
The wetted depth can be increased up to 56 cm (22 inches), The test section is 
bounded by glass windows on sides, bottom and top to permit visual and optical 
study of the flow throughout its 1.83 m (6 feet) length. 
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The existing arrangement made it possible to achieve a maximum velocity of 
1 m/sec in the test section. Turbulence levels and velocity profiles at 20 cm 
from the entrance to the test section and outside of the wall’s boundary layer 
were obtained using hot film probes and an LDV system. The uniformity of the 
flow within a square area (30 cm x 30 cm) centered at the center of the test sec- 
tion was below 1% of the center line mean velocity. The free stream turbulence 
level at 25 cm/sec on the centerline was near .05%. The power spectrum of 
free stream velocity fluctuations, a sample of which is shown in figure 2a, indi- 
cates most of the energy was concentrated at frequencies below .5 Hz, whichare 
primarily below the range of frequencies of interest in this study, 2 to 10 Hz 
(figure 2). 

2 . 2 . Model 

An axisymmetric model, constructed of anodized aluminum and Plexiglas, 
was used throughout the investigation. The model was composed of a rear por- 
tion and a front portion which were connected by a sting as shown schematically 
in figure 3. The rear portion, a hollow cylinder 10.16 cm (4 inches) in diameter 
and 91.4 cm (36 inches) in length, was connected by a strut to a traversing 
mechanism which was installed on top of the test section. This allowed for 
movement of the model up and down and also along the test section axis. The 
front portion consisted of a half portion 2:1 ellipsoidal cone with a 10.16 cm (4 
inch) minor diameter. The ellipsoidal cone was extended in the downstream 
direction by a straight section which also was 10.16 cm (4 inches) in diameter 
and 5.08 cm (2 inches) in length. 

The sting was made from a 1.27 cm (1/2 inch) O.D. thick wall stainless steel 
tube. Twelve brass tubes were installed inside of its walls to carry pressure sig- 
nals from the front nose while the interior of the sting was devoted to dye tubes 
and electrical wires for the strip heaters. The front end of the sting was fixed to 
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the inside wall of the front nose while the rear end was connected to a vertical 
rod by means of a 90° gearing mechanism. The vertical rod was then brought 
out of the model through the strut. By manual rotation of the vertical rod, the 
sting could be moved longitudinally relative to the rear portion of the model. 
With the aid of this mechanism, plus a dial and a helipot at the end of the verti- 
cal rod, it was possible to change and read the cavity length. The solid inner 
faces of the cavity were equipped with 14 pressure taps on the downstream face 
and 8 pressure taps on the upstream face. The diameter of each tap opening 
was .1 cm (.041 inch ). 

2.3. Pressure Measurements 

The low velocities in the range of interest introduced the well known difficulty 
of measuring low differential pressures under high head. A free stream velocity 
of 25 cm/sec generated a dynamic head of .38 cm (.15 inch) of water, while a 
typical value of total head at the point where the pressure transducer was 
located was 63.5 cm (25 inches) of water. Further, since differential pressures 
inside the cavity are much smaller than for a free stream dynamic head, the 
demand on the pressure measuring system for high resolution and repeatability 
is very severe. 

Another common problem in pressure measurements with water as a medium 
arises from the presence of air bubbles in the system, which may cause errone- 
ous pressure readings. The following approach was taken to solve the air bubble 
problem and to achieve high resolution. Each line, Tygon tube ID = .16 cm (1/16 
inch), from the model pressure taps or from the pitot-static probe was con- 
nected to a common drain through a toggle valve (figure 4). Each line then 
extended and manifolded into a 24 port fast switch (Scanivalve type W2460). 
The common end of the Scanivalve was connected to the low pressure side of a 
pressure transducer (Validyne type Wet-Wet DP-103-12) which was able to 
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measure differential pressures in the range of 0 to .762 cm (.3 inch) of water 
with a resolution of 1% of the total range. The high pressure side of the trans- 
ducer was connected to the total pressure tube in the test section, which was 
chosen as the reference. The common drain was connected to a vacuum line and 
was used as a sink into which air bubbles could rise. Also, it was used to fill the 
lines from the model to the toggle valves. A special calibration tank was 
designed to provide a differential head of 0 to 10.16 cm (4 inch) of water with a 
resolution of ,00254 cm (.001 inch) of water. One port of the Scanivalve was 
devoted to the low head side of the calibration tank, while the higher head side 
was connected to the high pressure side of the transducer through a 3-way 
valve. A schematic of the pressure sensing system is given in figure 4. The aver- 
age height of the calibration tank was set close to that of the free surface. This 
was done to prevent overpressurizing the transducer during switching. Table 1 
gives different operational positions of the valves. The pressure measurement 
process was automated by using an HP data logger. A proper software program 
was written for the data logger in order to actuate the Scanivalve and to digitize 
the analog signal from the transducer. Prior to each run the digitization rate, 
the lag time between each point, and the averaging time were specified to the 
machine. After each run, the program stored the pressure distributions, and 
the calculated drag coefficient as well as other parameters of each run on a 
floppy disk. 

2.4. Frequency. Phase and Correlation Measurement s 

Velocity fluctuation measurements were made with a TSI hot -film probe 
(1 277-1 0A) using a TSI 1050C anemometer in conjunction with a TSI true rms 
meter (1076). The output of the hot film anemometer was fed to a digital power 
spectrum analyzer (SL20). A bandwidth of .1 Hz was used for locating and dis- 
tinguishing frequency content of the velocity signals on line. A hard copy of the 
spectrum was obtained by either taking pictures of the spectrum analyzer 
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screen or by plotting the spectrum on an x-y plotter (HP model 700 1A). 

For phase angle and cross correlation measurements a second hot-film probe 
was mounted on a traversing mechanism which was able to traverse the entire 
width of the cavity. Throughout this experiment the first probe was used as a 
reference probe. This probe was placed above the upstream corner at x = b and 
y. a5 , the latter being the y-location of u/U e = .95. Phase measurements were 
made by moving the second probe relative to the reference probe along the edge 
of the shear layer where U/U e = .95 ( y = y % ss ). The signals from the reference 
and the second probe were fed to a digital probability and correlation analyzer 
(SAICOR 43A). Throughout the phase angle or cross correlation measurements a 
sampling rate of 200 per second, limited by 400 samples and a corresponding 
Nyquist frequency of about 30, was used. Again a hard copy of the output was 
obtained on an x-y plotter. See Appendix. A for the details of phase angle meas- 
urement. 

2.5. Velocity Measurements by Laser Doppler Velocimetry 

Due to sensitivity of the cavity flow oscillations to the presence of external 
objects in the flow, hot-film anemometry is not suitable for velocity measure- 
ments. Severe interference of the probes with cavity flow mode switching is 
observed in the present study. Existence of large regions of reversed flow is 
another reason for making hot-film anemometry unappealing for this flow. The 
only instrument currently available to measure velocity field non-intrusively and 
unambiguously in such a flow condition is the laser Doppler velocimeter (LDV). 
Other useful properties of the LDV include the small sampling volume, high fre- 
quency response and linearity of the instrument. In this study, the velocity field 
of the cavity flow was measured by a two-component, two-color, frequency 
shifted laser Doppler velocimeter (TSI 9100-7 system). It is a dual beam (for- 
ward) scatter configuration and processor, which was operated in single-burst 
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mode. The LDV processor operated in conjunction with a minicomputer system 
(PDP-11), which was also used for processing the results. A forward scatter 
configuration was used in these experiments, with the transmitting and receiv- 
ing optics located at the bottom and top side of the facility as shown in figure 5. 

The physical arrangement of the LDV optics relative to the facility and infor- 
mation regarding focal volume are given in figure 6 and table 2. The transmit- 
ting and receiving optics are mounted on an optical table with their optical axes 
perpendicular to the flow direction. The focal volume has its longitudinal axis in 
the z direction of the facility, The streamvri.se velocity component is measured 
by the fringe system formed in the z-x plane of the tunnel >while velocity com- 
ponents perpendicular to the mean free stream are measured by the fringe sys- 
tem formed in the z-y plane of the channel. The velocity bias technique (Bragg 
cell) is used to detect velocity reversals. For details of the optical setup and 
Bragg shifting see Appendices B and C, respectively. 

The focal volume can be moved within the test section by means of remotely 
operated stepping motors. Motion on a fixed horizontal plane was obtained by a 
corresponding motion of the optical table. Motion along the z-axis resulted 
from the displacement of the focussing lens and receiving optics along their opt- 
ical axis. The table position was controlled to a resolution better than 100 n m. 

2.5.1. Signal Processor and Data Acquisition. The measurement of the velo- 
city can be obtained only after suitable processing of the outputs of the pho- 
tomultiplier tubes. Two counter processors were used, one for each velocity 
component. They were operated in the single-burst mode, for which only one 
particle is present in the focal volume when the measurement is made. In addi- 
tion, the system was operated in coincident mode, in which the output of the 
processors is tested for coincidence in time so that the two components of the 
velocity vector as determined by the optical set-up are measured on the same 
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particle. 

The signal processing is accomplished in four stages: analog processing, burst 
processing, coincidence test and transfer to the computer. The analog process- 
ing involves suitable filtering and amplification of the signal for optimum signal- 
to-noise ratio. Only low pass filters were used. A built-in filter in the Bragg cell 
circuitry removes pedestal noise (low frequencies). The low pass filter is used 
for removal of high frequency noise. Cut-off values from 30 kHz to 100 kHz were 
used, depending on flow velocity. The output of the filter is processed by an 
amplifier. In the present experiments the gain was set at 2 or 4 depending on 
the condition of the signal. A detailed description of the valid burst processing 
is given in Appendix D. 


2.5.2. Operating Conditions. A typical signal-to-noise ratio, as observed from 
high sampling rate measurements in a quiescent fluid, was in excess of 100. An 
equivalent background noise level of .4% and .6% was detected for the stream- 
wise and the transverse velocity fluctuations respectively. When the signals 
from two channels were correlated for the purpose of shear stress measure- 


ment. the corresponding noise level for 



measurements was ,07%. 


The 


better noise level value was due to the fact that much of the uncorrelated noise 
from two independent processors was eliminated during the correlation process, 
resulting in an extremely low noise level for shear stress measurement. 


A time window of 50 psec was used in all velocity measurements. The ampli- 
tude threshold level and the gain of the amplifier were adjusted to maintain the 
data rate between 60 to 200 particles /sec, which corresponds to a minimum 
Nyquist frequency of 10. It was found under the present experimental condi- 
tions that particles less than 15 had sufficient response frequency. Thus Jor 
the purpose of LDV measurements, particles larger than 15 pm were filtered 
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from the water. 

2.5.3. Operating Procedures . Prior to every experimental run with the LDV 
system, a variety of tests were conducted to verify adequate operation of the 
system, e.g., the transmitting optics alignment and overlapping of the green and 
blue focal volumes. The alignment of the receiving optics was verified with the 
facility in operation. Adjustment of the location of the pin-hole and photomulti- 
plier assembly was often necessary for maximum data acquisition rate at the 
operating conditions of the water channel. 

After the high-pass and low-pass filters were set to their predetermined 
values, the gain of the amplifier stage was adjusted for optimum signal-to-noise 
ratio. After these tests were completed, the actual measurements were ini- 
tiated. A 16,000 word buffer was obtained at every point. The data were stored 
temporarily on disk. At the end of the test run, all the data were transferred to 
a magnetic tape for permanent storage and processing at a later time. The size 
of the data buffer combined with the data rate to give a total sampling time 
between 15 and 60 sec. This sampling time was considered adequate to properly 
characterize mean turbulent quantities. 

2.5.4. Data Processing. The velocity data generated by the data acquisition 
programs was processed using the minicomputer system. Processing involved 
several steps, e.g., error detection, computation of flow characteristics in the 
frame of reference of the optics, transformation of those quantities to the phy- 
sical frame of reference, and the computation of integral quantities in the phy- 
sical reference frame. 

There was one condition that could have led to a significant error in the 
measurements. This condition arises from transmission line noise, resulting in 
a spurious data word being included in the stream of data as it is transferred to 
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the computer. Because catastrophic effects could, result from this condition, 
the data processing programs had to verify that the data were in the correct 
order. If an error was detected, the data point would be discarded. Less than 
1% of the data was removed by this test. At every measurement point, a fre- 
quency histogram was constructed for each channel. This histogram provided 
information about the proper operation of the system, i.e., filter setting and the 
electronic noise. 

The Doppler frequency and the time between data points were utilized to cal- 
culate the flow velocity using the appropriate Doppler constant. The mean, 
second moment and cross-correlation between the two velocity components 
were evaluated. For a detailed description of the expressions used for the data 
processing see Appendix E. 

2.8. Flow Visualization 

The model was equipped with two dye injection tubes whose axes each made 
an angle of 1° with local curvature of the nose cone at the point of the injection. 
The injection points were 2.54 cm (l inch) off the stagnation point. Two visuali- 
zation techniques were employed in this study. In both cases, diluted coloring 
dye or fluorescein dye was gravity fed to the injection tubes and then discharged 
into the boundary layer over the cavity. 

In the first method, using blue dye, a uniform background illumination was 
provided specifically for photographic purposes. In the second technique 
fluorescein dye was used and excited by a laser sheet. The laser sheet was 
obtained from a 2-watt argon-ion laser with a thickness less than 1mm. The 
technique was used first by Dewey (1976). Photographs of the resulting flow pat- 
terns were made using a 35 mm SLR camera with a focal length F3.6 lens at 
shutter speed 2 to 1/1000 sec and either TEI-X black and white print film or 
high speed Ektachrome tungsten color slide film. 
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2.7. Operational Procedures. Accuracies and flow Parameters 

The pressure and velocity measurements procedure was as follows. Prior to 
each run with the tunnel running at idle, the nullity of the pressure readings 
was checked. The desired speed of the tunnel was set. The tunnel speed, 
extracted from the dynamic pressure measurement, was checked against the 
one measured by laser Doppler velocimeter. The cavity gap was set to the 
desired width and the frequency of cavity flow oscillation was checked. Depend- 
ing on the degree of the flow unsteadiness, the pressure averaging time was set 
from 30 to 120 seconds. The alignment of the model was checked by pressure 
reading differences of taps on the four sides of the downstream face. A cross 
check of the model alignment was made through flow visualization. 

At the end of a test run, with the tunnel again set at idle, the tap pressures 
were checked to insure they had returned to zero. For velocity field measure- 
ments, the measuring volume was fixed at a reference point on the model. This 
was done by observing the position of measuring volume with respect to the 
model through the eyepiece on the receiving side of the LDV. Measurements 
were made at speeds covering the range of 16 to 27 cm/sec. The corresponding 
Reynolds number, Re, based on the model diameter, d = 10.16 cm (4 inches), and 
at U B = 23 cmAec was 24 x 10~ 3 . Most of the pressure data points were meas- 
ured at least twice and at some points up to 5 times. The repeatability of the 
calculated drag coefficient was about 4% of the average value. 

2.B. Reference Quantities 

Throughout the experiment, the origin of the reference coordinate system 
was located at the upstream corner, The x and y axes lie in the symmetry 
plane. The x-axis was in the direction of flow; the y-axis was measured positive 
outward with respect to the center of the cavity. The components of velocity 
vector in the direction of the x and y axes were denoted u and v. 



- 1B- 


The cavity width, b, was non-dimensionalized by normalizing it with respect 
to 6q , the momentum thickness of the upstream boundary layer at the point of 
separation. Momentum thickness for the shear layer is defined as 


4-w 

e=f — 

J u B 


'-■h dy 


( 2 - 1 ) 


where C/„ is the free stream velocity on the outside edge of the boundary layer 
at the point of separation. Using U g instead of [/„, the free stream velocity in 
front of the model nose eliminated the need for the blockage correction. This 
also made the results independent of the front nose shape (provided there 
existed no substantial pressure gradient over the length of the cavity). 

The boundary layer upstream of the cavity was considered laminar at 

6 * 

U g = 23 cm /sec and had a shape factor of H = — = 2.5 where 5* = .1 cm (dis- 
placement thickness), fib = 0.041cm and Re^ = 95. For the shear layer, the 

[y - yJ 

transverse distance y was non-dimensionalized as , where y, s is the 

location of the non-dimensional velocity — — with a value of .5, and e is the 

local momentum thickness of the shear layer as defined by equation (3.1). The 
radius of the axisymmetric cavity, which is 5,08 cm (2 inches) was used to nor- 
malize the radial distance of pressure taps (r) on the inner surfaces of the cav- 
ity. All drag forces were normalized by ^-pU g A , where p is water density and 
A is the cross-sectional area of the model. 
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Chapter 3 

Flow Documentation 

The measurements of various parameters relevant to the cavity oscillations 
reported in this chapter are not all novel, some of them having been reported in 
references cited above and in others. The purpose here was to provide a solid 
reference condition for the measurements to be described in Chapters 4 and 5. 
Nevertheless, some new definitions, clarifications and features are intro- 
duced . 

3.1. Cavity Plow Oscillations as a Function of b/e 0 

It is known [Karamcheti (1956), Sarohia (1975), Rockwell and Naudascher 
(1979)] that shear layers over cavities undergo transition from a steady layer to 
an unsteady layer at a critical dimensionless width b /6b I min • For 
b/e o < b/9o I min. the shear layer bridges smoothly over the cavity with no distinct 
oscillation in it. The onset of unsteadiness is always associated with the presence 
of a distinct fundamental frequency, which is amplified by the shear layer. In 
this study, the fundamental frequency is identified as the one that has the 
highest peak in the power spectrum of the velocity fluctuations. Velocity 
fluctuations were measured at x = b and y =v. 95 .(v .95 being the locus of 

yj—= .95) throughout the frequency measurements. 

Figure 7 presents the variation of fundamental frequency, F, with b/% at a 
fixed edge velocity {U e = 23cm /&ec). Oscillation starts at b/ 0 o « 81, which 
represents the minimum cavity width 6 / 00 !,^, for the present experimental 
condition. The fundamental frequency of oscillation decreases as b/e 0 
increases. This course of frequency variation is interrupted by a sudden jump 
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of frequency from F = 4.6 to 6.2 Hz. This new frequency initiates another 
course of decreasing frequency as 6/% is further increased. A non-dimensional 

frequency can be defined for these two modes of oscillation as — — (Strouhal 

Number). This Strouhal Number has a value of « 1 over the entire range of the 
starting mode, then, jumps to a value of « 1.5 which it maintains over the next 
mode. According to Sarohia (1975) and Rockwell and Knisely (1980) these two 
modes correspond to Modes II and III of the cavity flow oscillation. 

A sequence of power spectra of the velocity fluctuations as a function of b /9 q 
is shown in figure 8. Several features of these spectra are noteworthy. First is 
the presence of peaks both at 4.6 Hz and 6.2 Hz during the transition region 
between the modes. From these average power spectra it is not possible to 
determine whether the two frequencies alternate intermittently or whether both 
coexist simultaneously. To clarify this matter, figure 9 shows three oscilloscope 
traces of the velocity fluctuations in Mode II, in the transitional region, and in 
Mode III of the oscillations, respectively. 

Figure 9b confirms the coexistence of two frequencies by showing a strong 
modulated trace of the signal. A region of hysteresis, as shown in figure 7, is 
associated with the transition region, Presence of frequencies lower than funda- 
mental is a prominent feature of Mode III (figure 8c). Such low frequencies have 
been observed by Knisely (1980) and were shown to be associated with the 
interaction of successive shear layer vortices with the downstream corner. A 
possible explanation for the existence of such a low frequency is given in the 
flow model section (Section 4.10). End of Mode III is marked by a loss of coher- 
ence in the oscillation and the disappearance of the fundamental frequency. At 
this stage (6/9 > 160) , cavity flow begins to exhibit bluff body wake behavior, 
and a distinct frequency for this regime cannot be detected. 
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3.2. Phase Measurements 

Information about the overall phase difference between fluctuations at the 
upstream corner (x = 0) and the downstream comer (x = b) is needed for 
improving our understanding of the cavity flow oscillation. This is the subject of 
the following discussion. 

The existence of a strong phase gradient across laminar and turbulent shear 
layers is well known [Freymuth (1966), Sarohia (1975), Elder (1978), Ziada and 
Rockwell (1981), and Hussain (1981)]. This gradient can be well approximated by 
inviscid theory (Michalke, 1965), except for regions close to the downstream 
corner. Existence of such a phase gradient introduces difficulties in choosing a 
proper criterion to define an overall phase difference between the two corners of 
the cavity. Kmsely (1980), Ziada and Rockwell (1961), Rockwell and Schachen- 
mann (198 2), in their effort to define an overall phase difference for different 
types of self-sustaining flows, concluded that tracking values of phase along the 
edge of the shear layer tu/U B = .95) will result in an overall phase difference of 
~ 2N ix . In contrast, other criteria.e.g., following one of the disturbance ampli- 
tude maxima within the shear layer, yield different overall phase differences for 
different geometries of self-sustaining flows. 

In the present study, phase measurements were obtained by following the 
outer edge of the shear layer along the locus of u/U e = .95. The reference probe 
was located at x = b, y = y t9s , while the second probe was moved along t/ > 95 of 
the shear layer. Before each set of measurements was to be taken, the auto- 
correlation of the reference signal with respect to itself or to that of the second 
probe (cross-correlation) was obtained. In this way, any initial phase difference 
due to circumferential or longitudinal distances between two probes was 
detected. A sample of such a measurement is given in figure 10. For the case of 
b/e 0 = 85 and U e = 23cm./&ec a typical phase angle distribution of the velocity 
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fluctuations as a function of x/e^ is presented in figure 11. The cavity shear 

Fb 

layer has a fundamental frequency of 5.9 Hz with corresponding — -= 1» which 

Uq 

identifies it as the second mode of oscillations. Figure 11 depicts an overall 
phase difference of 4 rr between two corners of the cavity. Similar measurements 
to determine the overall phase difference between the two corners were carried 
out for the entire range of Modes II and III of cavity flow oscillation. Figure 12 
represents the distribution of the overall phase difference as a function of b/9 0 . 
It is important to note that the cavity flow retains a 4 tt phase difference during 
Mode II, while having a 6 tt phase difference during Mode III. This fact suggests 
that overall phase difference between the two corners be formulated as 


^-=N 

2 7T 


(3.1) 


where N represents the mode of the cavity flow oscillation. Using the relation 
U c = F\ , where U c is the convection velocity (phase velocity) of the fundamen- 
tal frequency of oscillation and \ is the wavelength, equation 3.1 can be rear- 
ranged as 


y _ b _ Fb 
2n ~ X = U c 


(3.2) 


Equation 3.2 indicates that cavity during Mode II contains two wavelengths of 
the fundamental frequency, while Mode III contains three wavelengths. Also, 
emergence of the Strouhal Number, based on the phase velocity in equation 3.2, 

Fb 

directs us to define the mode of the oscillations based on — - = N rather than 

k'c 

Fb N 

___ = __ Thi. s convention defines each mode by the number of wavelengths 
Uq s 


that it contains and establishes a relation between the mode's name and the 
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corresponding Strouhal Number. 

It Is important to note that close to the transition region, an important phase 
deviation from « 2 occurs (figure 12). Such a deviation of overall phase 

difference also was observed by Rockwell and Schachenmann (1982) for imping- 
ing turbulent jets at low Mach numbers. 

3.3. Overall Features of Cavity Flow Regimes 

3.3.1. Nan-oscillating Cavity Flaw. The non-oscillating regime of the cavity 
flow is marked by a fairly steady and weak vortex stationed inside of the cavity. 
Unsteady mass transfer in and out of the cavity is minimal due to lack of any 
substantial oscillations in the shear layer. Figure 13 presents flow visualization 
of the cavity flow in this regime. Often, color dye trapped in the cavity remained 
there for a long time. The effect of the cavity presence on the outer flow field is 
minimal. Branching of the shear layer at the downstream corner is another 
feature of this mode observed in flow visualization studies. It should be men- 
tioned that for b/e 0 < b/8 0 | ^.hot-film signals indicated the presence of weak 
periodic fluctuations in the shear layer. Due to lack of a wide enough cavity , 
these perturbations do not amplify enough to make the flow self-sustained. This 
matter will be discussed in more detail in Chapter 5. 

3.3.2. Self-Sustaining Cavity Flow. The definition of self-sustaining flow and 
the description of the chain of events that result in flow self-sustainment is 
given in the introduction. As mentioned before, natural cavity flow oscillation is 
by definition self-sustaining. 
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3.3.2. 1 General Characteristics of the Self-sustaining Modes.. A self- 
sustaining oscillation appears in the cavity when b/8 0 reaches b/Soimin- The 
wavelength of the fundamental frequency scales itself with the width of the cav- 
ity (figure 12). Plow in this regime is stable. Even severe attempts to perturb 
the course of oscillations is resisted, depending on the strength of the self- 
sustainment. In this study, the most stable case is defined as the one that gen- 
erates the highest level of organized fluctuations for a given b/e Q . 

Overall phase difference between the two corners of the cavity satisfies the 

relation -^— = N, where N is the number of wavelengths contained 

2rr A U c 6 

by the cavity width during Nth mode of oscillations. Shown in figure 14 are flow 

visualizations of Modes II and III of the cavity flow oscillations. In the present 

study, the cavity always started the oscillation with Mode II, omitting the first 

mode — — ra 1. Sarohia (1975) observed similar omission of the first mode in 

u c 

his study of cavity flow for some free stream velocities. 

Instantaneous flow visualization pictures (figure 14) show that the wave 
amplitude grows as it approaches the downstream corner. This growth, which is 
exponential according to the linear theory, results in a vortex roll up. This roll 
up is very similar to that observed by Brown and Roshko (1974), Winant and 
Browand (1974) in the free shear layer. Mode III shows, usually, two complete 
roll ups of the shear layer. Unlike the case in a free shear layer, no vortex pair- 
ing is observed during the cavity self-sustaining modes. Thus, presence of low 
frequency peaks in the power spectra of Mode III of the present study cannot be 
associated with vortex coalescence. 

Presence of a strong stable vortex is another characteristic of the self- 
sustaining mode. 
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3. 3.2. 2 Similarity to Externally Forced Free Shear Layers.. The ladder-like 
variation of cavity shear layer frequency shows striking resemblance to the 
externally forced free shear layers (Miksad 1973, Ho and Huang 1982). One can 
imagine that the downstream corner plays the role of external forcing source. 
Also the term bAJ c can be assumed to be equivalent to the forcing frequency in 
an externally forced free shear layer experiment. Figure 15 shows the varia- 
tions of the fundamental cavity shear layer frequency as a function of b/U c . In 
this figure, several cases of the cavity flow geometries and velocities fall on the 
corresponding mode lines of externally forced experiments. In general, one 
must treat this resemblance with caution, because at high subsonic or super- 
sonic cavity flows, the wavelength will not scale with cavity width b. As sug- 
gested by Laufer and Monkowitz (1980) such a scaling should follow the more 
complicated feedback equation 


Fb 


+ 



(3.3) 


where C is the sound velocity. In the subsonic case C» U c , equation (3.3) 

Fb 

reduces to our familiar relation for the Strouhal Number ——• = N. 

u c 

3.3.3. Wake Mode of Cavity Flaw. The last Mode (III) of oscillations, observed 
in this study, shows strong intermittency as the cavity width is increased. It 
seems that eventually an intermittency overcomes the coherent oscillation and 
the feedback mechanism becomes ineffective. The flow visualization studies 
revealed that when the self-sustaining regime ends, the shear layer no longer 
reattaches onto the downstream corner. The cavity flow becomes semi- 
periodically unstable in a larger scale, comparable to the scale of the model 
(figure 16a, b). Time -averaged flow visualization (figure 16c) shows that the 
mean flow reattaches below the downstream corner. Instantaneous flow 
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visualization indicates that the flow is highly asymmetrical (figure 16a). The 
flow field m this re gim e shows some resemblance to a three-dimensional bluff 
body wake (Werie, 1974), even though the presence of the rear body is still 
strongly being felt by the wake of the front body. A detailed exploration of the 
flow for cavity-diameters ratio b/R > 1.65 was not carried out. 


3.4. Mean Velocity Distributions in the Shear Layer 

For three values of cavity width b/6$ = 65 ,b/6o = 103 and b/ 0 Q = 130. mean 
velocity profiles were obtained at a number of streamwise stations. Edge velo- 
city (U e ) was fixed at 23 cm/sec. The cavity width b/e^ = 85 marks the begin- 


ning region of Mode II 


Fb 

U c 


= 2 


middle of the same mode. For b/6 0 = 130 


, while the cavity width b/e^ - 103 marks the 
Fb 


U c 


-= 3 


. the cavity was oscillating in 


Mode III. Figure 17 presents nondimensionalized mean velocity profiles at 
selected downstream stations for 6 /fife = 05. As shown in figure 17. in the early 
stages of the shear layer growth, the velocity profile changes from a boundary 

\L 

layer profile to a shear layer profile. Free-shear-layer similarity m the — — 

Uq 

profiles is attained for b / > 30. In figures 18 and 19, which present non- 
dimensionalized mean velocity profiles for b/e o = 103and&/6b = 130, respec- 
tively, self-similarity in the yy— profiles is also observed. 


Figures 20 through 22 present the nondimensionalized mean transverse com- 
ponent of the velocity, y-—, in the shear layer. These figures reveal that unlike 

Ue 

the -jy— profiles, mean transverse velocity profiles do not scale with the down- 
Ue 

stream distance. Thus>no self -similarity can be found for yr- profiles. Gorree, 

Ue 

Dimotakis and Koochesfaham (1900) were able to show that for free shear 
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V 

layers with very high Reynolds numbers, self-similarity exists for profiles. 

Uq 

"V 

The lack of similarity in the --—profiles of the cavity shear layers might be due 

Uq 

to existance of the strong recirculating flow inside of the cavity. Figure 23 

X 

presents the locus of y $ values as a function of — As in previous studies by 

Browand (1966), Sarohia (1975), and Knisely (1980), the y, 5 values were found 
initially to lie above the level of the separation edge. Depending on the stage of 
the cavity flow, the locus of y . 5 approaches the y-0 line as the profile 
develops. The shear layer shows a stronger deflection into the cavity for the 
b/e 0 = 103 than for b/e = 85. 

3.5. Growth Rate of the Shear Layer 

For four values of cavity width, b/6o = 85 (Mode II), fe = 103 (Mode II), 
b/0Q = 119 (Mode II) and b /6 q = 135 (Mode III), the momentum growth of the 
shear layer is given in figure 24. As in the studies by Sarohia (1975) and Knisely 
(1980), the momentum thickness is found to grow linearly with x/e^. Because of 
this linear growth, one might assume a self-similarity for the cavity shear layers, 
but it is interesting to mention that in both the previous studies and the 
present experiment, the value of the shear stress is not constant (figure 24) as 
is required for the self -similarity of the free shear layer. 

d, 0 

For all cases with 6 /©o ^ 103, the growth rate is fairly constant ~^T = «031 

and does not experience a sudden change during the mode switching region. It 
should be noted that the cavity shear layer growth rate is comparable to that of 

Ui 

turbulent mixing layers. For - 7 — = 0, the high Reynolds number shear layer of 

Uz 

Liepmann and Laufer (1947) had a growth rate momentum thickness of 0.035. 
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3.6. Shear Stress and Velocity Fluctuation Distribution in the Shear Layer 

Measurements of shear stress and fluctuating components of the velocity 
field are essential to an understanding of the overall nature of the cavity shear 
layer. As discussed previously, due to the nature of the present LDV measure- 
ments, shear stress measurements show fairly low background noise levels. 

Such an advantage was not available for -y-— or ~n~ measuremen t s - (figure 25) 

Figures 26 through 28 depict the streamwise evolution of the transverse 
profile of the shear stress for b/^) = 85 (Mode 11), b/6fe = 103 (Mode 11), and 

X 

b/e - 130 (III Mode). For large — ; shear stress profiles tend to show a second 

So 

peak in the lower portion of the shear layer. Existence of such a second peak 
has been observed [Freymuth (1966) ,Knisely (1980) ] ; Ziada and 'Rockwell 
(1981) showed that it can be predicted by Stuart's (1967) theoretical vortex 
model. Also important for our consideration is the fact that, even though the 

X 

mean streamwise velocity profiles scale with — -, the maxima of the shear stress 

e o 

profiles do not have a constant value as would be required for self-similarity, 
they show an exponential growth with x/Sq (Figures 29 and 30). Streamwise 
growth of the maximum shear stress for several cavity widths is shown in figures 
29 and 30. Short cavity cases, b/do = 85 (Mode II) and b/6a = 103 (Mode II). 
show a region of exponential growth of the maximum shear stress. For long cav- 
ities b/&Q > 103, the region of exponential growth is followed by a saturation 

region. In Mode II, the rate of growth for — g — increases as b/&o increases 

, U e 

(figure 29). The transition region does not show a noticeable change in the 
growth rate, but the general level of the shear stress decreases for the transi- 
tion region (figure 29). 
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Figure 30 gives the distribution of maximum shear stress growth for Mode 111. 
This mode, too, shows exponential growth for the maximum shear stress m the 
shear layer. Agam,growth rate increases as b/&o increases. It is important to 
mention that the general level of the maximum shear stress distribution is 
higher for b/e^ = 151.8 and b/6$ = 158.5 .which are close to the end of Mode III. 

The maximum shear stress achievable during Mode II, 0.0086, is lower than 
that in Mode III. which is .013. For comparison the value of maximum shear 
stress for a self-similar, two-dimensional turbulent mixing layer as determined 
by integrating the velocity profiles presented in Liepmann and Laufer (1947) is 
.0125. Integration of the velocity profiles measured by Wygnanski and Fiedler 
(1970) gives r majL /pU z e =. 015. Kistler and Tan (1967) reported a maximum 
value of 0.012 for two-dimensional cavity shear layers. 

Figures 31 and 32 depict the variation of the maximum streamwise and 
transverse fluctuation as a function of 6 /% Again an exponential growth of the 
disturbances is evident. It is clear from these figures that rms fluctuations 
suffer from a high percentage noise level of either flow background or electron- 
ics. 


3.7. Mean and fluctuating Velocity Field of the Cavity Flow 

Figures 33 through 37 present traverses of the velocity field through the cav- 

’ Fb 


ity for one selected cavity width ( b/e^ - 103) m Mode II 


V. 


-= 2 


The particu- 


lar case selected has a fundamental frequency of 5.1 Hz. Figure 33 presents the 


mean streamwise 


u 


U B 


velocity profiles of the upstream boundary layer at the 


upstream corner, at four stations within the cavity and two stations down- 


stream of the downstream corner. Similarly figure 34 displays the mean 


v 

U B 


profiles. 
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The change in flow direction inside the cavity in figures 33 and 34 is due 
to the recirculating vortex flow visible m the flow pictures. The upstream 
boundary layer is laminar and has a shape factor of 2.5. The downstream boun- 
dary layer profiles at x/d a = 110 and x/6 0 = 142 give respective boundary layer 
shape factors of 1.54 and 1,46. These shape factors are close to that of a tur- 
bulent boundary layer. 

Figures 35 and 36 reveal that significant fluctuations occur inside the cavity, 
but the flow field inside the cavity has a much lower shear stress level than that 
of the shear layer (figure 37), possibly revealing, that most of the u' and v' 
fluctuations are not well correlated. Perhaps the most important difference 

' u ' 

between — — and --— behavior is shown in the downstream boundary layer 

profiles (figure 38). The maximum value of the yj— profiles in the shear layer 
continues to increase, even after the shear layer transforms to a boundary 

'XL 

layer. On the contrary >the maximum value of the -jj— profiles shows a decrease 

after the downstream corner. The shear stress profile shows a monotomc 
decrease in the maximum value of the profile as the shear layer transforms to 
the boundary layer. 
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Chapter 4 

The Drag Coefficient of the Cavity 

Throughout this study the term "drag coefficient of the cavity" refers to the 
net force in the flow direction being experienced by the gap. The overall changes 
in the net drag coefficient of the model,which might arise from cavity presence, 
will not be discussed in this study. 

4.1. Momentum Balance of the Cavity 

A momentum balance for the fluid enclosed by the control surface shown in 
sketch 2 gives 

/ JLUL’ nJdAi = -/ P mlAz A + frn.dAi S (4.1) 

s s $ 

where dAi = nddx^ , dAz A — Znrdr, dAg = nddx% and SZ_ = ui_+vj_, The x- 



SKETCH 2 
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component of momentum balance will be 

+ pnd f uv dx 1 = -2jt f - P^rdr + ndf T xy dx 1 + T>d a f T^dx^ (4.2) 

Introducing u = u + u' ,v = v + v' to eq. (4.2), taking the time average and 
neglecting the viscous shear stress on the free surface will give 

+prr d f tjuv + •ZfiT)dx 1 = -2 n f \p z - P 4 Jrdr + 7rd f dx 3 (4.3) 
Using the following definition will reduce equation (4.3) to equation (4.4) 


Cjjp = cavity pressure drag coefficient = 2ir 


f[p 2 ~ Pa rdr/~~p f/ e 2 rrd 2 /4 


C DS = drag coefficient of the sting = —pnd s Jr yx dx 3 /^pU a rrd 2 /4 


—pnd ftuv + u r iT)dx 1 = + Cpp + Cps 


(4.4) 


However, by definition, the drag coefficient of the cavity can be written as 

Cp - Cpp + Qw (4.5) 


The combination of equations (4.4) and (4.5) results in a new relation for the 
drag coefficient of the cavity: 
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Cj) — — rrd f (©17 + ©V)da: 1 /^-/yu z 8 rrc£ 2 /4 


(4.6) 


Relation (4.6) reveals that it is possible to determine the cavity drag coefficient 
by evaluating the mean and turbulent momentum transfer terms along a 
straight line connecting two corners of the cavity. Throughout this section, this 
line will be called y = 0 line, Roshko (1955) found that the skin friction 
coefficient along the walls of two-dimensional cavities were two orders of magni- 
tude less than the pressure and pressure drag coefficient. Therefore, skin fric- 
tion in the cavity as a contribution to the drag is entirely negligible and 

C D * C DP (4.7) 


4.2. Measurement of the Cavity Drag Coefficient 

The drag coefficient of the cavity is determined by measuring the integral of 
the pressure distribution over the solid walls of the cavity. Figure 39 presents 
the variation of drag coefficient with b/9 0 in. the three distinct regimes of cavity 
flow as defined in chapter 3. This figure reveals two rather significant charac- 
teristics of the cavity drag. First, the average level of the drag coefficient stays 
low as long as the flow remains in either the non-oscillating mode or the self- 
sustained oscillating mode. A second feature is the existence of a critical b/do, 
in which, by exceeding its value, the drag coefficient of the cavity jumps to a 
value which is one order of magnitude higher than the typical average drag 
coefficient during oscillating modes. During the course of the wake mode the 
drag coefficient increases monotonically and reaches a value of 0.39 at 
b/e 0 = 250. No more measurements were made beyond this point due to model 


limitations. 
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4.2.1. Cavity Drag Coefficient in Non-Oscillating Mode. The drag coefficient 
in non-oscillating mode has a typical value of 1,3 x 10 -4 at b/e 0 = 77. When the 
cavity was replaced by a solid surface, the drag coefficient of this surface based 
on the upstream boundary layer shape factor was estimated to be .7 x 10 -3 . In 
turn, this drag coefficient is equivalent to that of b/e 0 = 98 of Mode II of the 
oscillations. However, based on observations in Chapter 3, the presence of the 
cavity reduces the shape factor of the downstream boundary layer and conse- 
quently may increase the skin friction drag on the rest of the model. At this 
point, it is not clear that a net drag reduction can be obtained for the whole 
model. 


4.2.2. Exponential Growth of the Cavity Drag Coefficient During Oscillating 
Modes. Distributions of the drag coefficient in Modes II and III, plotted on a log- 
linear scale, are presented in figure 40. At the beginning of the oscillating 
modes (Mode II in this study), the drag coefficient has a value of 0.9 x 10~ 3 at 
b/e 0 = 85. The drag coefficient then increases exponentially and reaches a peak 
value of 0.02 at b/e D = 108. Toward the end of Mode II, the drag coefficient 
decreases and reaches a minimum value of 0.004 at b/9 0 = 120. The transition 
region between Mode II and Mode III is marked by an increase in the drag 
coefficient value; C 0P = .0072 at b/G 0 = 128. Mode III exhibits another course 
of exponential growth for the drag coefficient value. A final value of 0.39 for the 
drag coefficient is reached at b/e^ - 250 before transition to the wake 
mode occurs. 


4.3. Cavity Pressure Distribution 

A better understanding of the drag generation mechanism can be obtained 
from pressure distributions along the solid surfaces of the cavity. 


Cp 


p zh 


Figures 41 through 43 present the pressure coefficient, 


, along 
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a radius of the upstream and downstream faces of the cavity for several 
selected values of b/e 0 . With each distribution the corresponding value of C D P 
is shown. 


4.3.1. Pressure Distribution for Non-oscillating Mode b/e^ < 82 . Figure 41a 
reveals that in the non-oscillating region, the cavity is slightly under pressure, 
i.e., C p < 0. No significant differences were observed between the upstream and 
the downstream face pressure distribution .which corresponds to closed and 
open symbols, respectively. 


4.3.2. Pressure Distribution for Oscillating Modes. 


1) Mode II 82 «s b/% < 120 


Figures 41b through 41 g show the cavity pressure coefficient distribution for 
Mode II of the oscillation. All cases have negative pressure distribution over the 
entire upstream face. The most important point to note in all cases is the steep 
pressure rise near the outer edge of the downstream face. In each case, C p 
reaches its maximum positive value (Cp mM ) at the corner. The magnitude of 
C D reaches a maximum at b/9 0 = 108, then decreases toward the end of Mode 
II. The downstream face experiences a negative pressure below the downstream 
corner with the same level and trend as that of the upstream face. Thus, the 
contribution to the total drag from these two portions of the cavity walls cancel 
each other. This leaves the outer edge of the downstream face with Cp > 0 , the 
main contributor to the cavity drag. It is also interesting to mention that in all 
cases, the average pressure on the upstream face is negative; and thus this face 
experiences thrust during Mode II of the oscillations. 


The pressure buildup near the upstream corner is evidently due to the 
deflection of the shear layer into the cavity. This deflection of the shear layer 
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causes a strong inflow to the cavity and contributes to the generation of a stable 
vortex. A large portion of the cavity experiences negative pressure coefficient 
due to the presence of strong recirculating flow inside of the cavity. Similar 
types of observations were made by Johannesen(1955) , Roshko(1955), McGregor 
and White (1970), Charwat, et al. (1961), Fox (1964) and Rossiter (1964). 

2) Mode III 120 < b/6 0 <£ 155 

Figures 42a through 42d present pressure distribution for Mode III. Except 
for a smoother pressure rise near the outer edge, the pressure distribution on 
the upstream face shows behavior similar to that of Mode II. 

4.3.3. Pressure Distribution for the Wake Mode. Pressure distributions for 
three cases of the wake mode are shown in figure 43 . Two main features dis- 
tinguish the wake mode from the oscillating mode. First, the entire upstream 
face of the cavity experiences positive pressure, i.e., C p > 0. Second, the stagna- 
tion point, the position of C Pmsg , steadily migrates toward the center of the cav- 
ity. We expect that for a large enough value of b, the stagnation point should 
settle at the center of downstream face. While pressure on the rear face is 
increasing, the average base pressure level decreases as b/9 0 increases (figure 

43 ). 

4.3.4. Overall Characteristics of the Cavity Pressure Distribution. Figure 44 
presents a collection of typical upstream corner pressure distributions for 
different regimes of the cavity flow. The start of the self-sustained oscillations is 
accompanied by an overall drop of the upstream face pressure distribution. The 
transition from Mode II to Mode III (6 /9 0 = 126) is marked by a substantial 
increase in upstream face pressure. A flow visualization study revealed 
that the absence of a strong recirculating flow inside the cavity might be respon- 
sible for the pressure rise of the upstream face. In Mode III, low upstream face 
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pressure recovers and shows a level similar to that of Mode II. 

Upon the loss of the self -sustained oscillations, another drastic decrease in 
the -average level of the upstream face pressure distribution occurs. At 
fa/S^ ^ giO (b/R « 1.7), the average upstream face pressure coefficient has a 
value of -0.13. 

Figure 45 represents variation of maximum pressure, Q, , with b/e 0 . The 
value of C 0 increases as 6 /Sn increases until it reaches a local maximum 

r max 

at b/eo = 108. This maximum is immediately followed by a valley which marks 
the last point of Mode II. The transition region and Mode III show a monotonic 
increase in the value of C p Appearance of the wake mode is marked by a sud- 
den increase in the value of CL from +0.17 to +0.43 at b/6o » 162, The value 
of C Pmai shows slight decrease during the wake mode. The postion of the stagna- 
tion point shifts toward the center of the cavity as b/e 0 increases. It is impor- 
tant that the location of C Pmia is fixed in the vicinity of the upstream corner in 
the self-sustained modes. Migration of this point toward the center of the cavity 
starts only when the cavity loses the self -sustained oscillation to the wake mode. 
The rather dramatic change in the pressure field of the cavity caused by exceed- 
ing a critical value of b/$o has been observed by Roshko (1955), Charwat (1961), 
i Fox (1968) and Maul and East (1963). However, the role of self-sustained 
oscillations in determining this critical b/6 Q was in general unnoticed. 

Variation of Cpp at the middle of two-dimensional cavities studied by Roshko 
(1955) and Fox (1968) as well as those of the present study are shown in 
fig. 46. It is interesting that in all three cases the location of steep 
jump in the value of agrees fairly well. 
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4.4. Observation of Sudden Drag and Pressure Change in Other Flows 

Several other investigators have reported a sudden change from a steady low- 
drag flow to a highly non-steady high-drag flow over cavities, cutouts, or bluff 
bodies in tandem. The intention in this section is to compare some of their 
results with those of the present study, and to give a possible explanation for 
the existence of the critical length, b/9 0 . 

Koenig (1978) in his study of flow over the combination of a disk and a 
cylinder, defined two regimes based on the drag behavior of the system, namely, 
sub- and super-critical. Drag in a subcritical regime has values one order of 
magnitude lower than that of a supercritical regime. Separated shear layer in a 
subcritical regime smoothly joins onto the rear body, while in the supercritical 
regime the shear layer becomes unsteady on a much greater scale. Koenig’s flow 
visualization study revealed that the optimum low drag case appears to be a 
more general example of flow over cavities. Figure 47 presents optimum drag 
measurements (Koenig, 1978) superimposed on the present drag measurements. 
Despite the basic differences in the geometry and the curvature of the shear 
layer, figure 47 confirms Koenig’s flow visualization conclusions. 

Koenig also reports the presence of a steady and stable flow field in the 
supercritical region. Production of sharp acoustic tones was another feature of 
subcritical regime that was reported by Koenig. The symptoms of subcritical 
regime are similar to those of simple>self-sustained>cavity flow oscillations. 

Mair (1965) reported a sudden change from steady low-drag flow to nonsteady 
high-drag flow near the base of a streamwise oriented cylinder as the gap 
between the cylinder and a downstream disk was changed. Zdravkovich (1977) 
observed an abrupt start of vortex shedding from an upstream cylinder accom- 
panied by an increase in drag of a system comprised of two cylinders in tandem 
when spacing between them was changed. 
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The presence of a feedback loop actively doing the task of cross- 
communication between two corners in the gap generated by two bluff bodies in 
tandem may possibly explain extra-stability and low-drag characteristics of such 
flow systems. 


The lock-on of the separated shear layer, or the position of Cp to the 
downstream corner due to feedback mechanism, has been overlooked by previ- 
ous investigators and has often been considered a secondary effect. A further 
discussion of this matter will be given in section 4.10. 


4.5. Estimate of C D by Evaluation of 


uv 

U* 


and 


err 

ui 


Terms Along y = 0 line 


An estimate of C D can be obtained by integrating the shear stress and mean 
momentum term over the periphery of the cavity (equation 4.6). LDV technique 
was used to measure the x and y component of the velocity vector, to determine 

-^-and terms along y = 0 line. For each case between 20 to 30, measure- 
U* Ug 


ment points were selected along the y = 0 line, depending on cavity length. 
Smaller spacing between measurement points was selected close to the down- 
stream corner. Due to the flow field symmetry of the wake mode, no repeatable 

*ufO u *u 

measurements could be done in that regime. Thus, evaluation of — 5- and — =- 

Uf Uf 

terms were limited to the oscillating regime of the cavity flow. 


4.5.1. Confirmation of Exponential Orowth for Cj). Figure 40 presents the 

7 7*9 f 1 fjrf 

variation of C D based on the integration of — a ~ and — g- along the y = 0 line 

ui uf 

with b/e 0. It is important to note that the exponential growth of Co during oscil- 
lating modes has been confirmed by this independent method. Except for the 
cases close to the end of Mode II and close to the transition region between 
Modes II and III, the absolute values of Co p and C D agree fairly well. Study of 
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the — — , —5- and —5- distributions along the y = 0 line is essential in 

U e U e U‘ U? 

further understanding of the mechanism which is responsible for the Cq 
behavior. This is the subject of the following chapters. 

4.6. Distribution of —— Along y = 0 lane 

Provided that the distribution of v at a fixed x position is uniform on the 
periphery of the cavity, the integral of v along the y = 0 line should be minimal 
for a satisfactory mass balance of the cavity. 

Starting from figures 48a through 48g, a sequence of distributions for 

Ug 

a few values of b/e 0 in Mode II is given. Each figure is accompanied by its 
corresponding drag coefficient. 

V 

A common feature of all —— distributions is the existence of a narrow region 

of inflow to the cavity H7 close to the downstream corner. The negative v 
region is balanced by a broader region of the outflow from the cavity +17 close to 
the upstream corner. The stronger the inshooting of flow into the cavity, the 
higher is the peak of the positive v region. The peak of -17 increases as b/ 0 o 
increases and reaches a maximum at b/e= 112.5, which corresponds to the 
maximum drag case. This peak eventually decreases toward the end of Mode II, 
In Mode III the v region becomes slightly broader and shows a small change 
(figure 49). 

In terms of stable vortex inside the cavity, distribution can be inter- 

k'e 

preted as the induced velocity of cavity vortex on the y = 0 line. Therefore, it 
can be concluded that the strength of the vortex increases during the Mode II as 
6/©o increases. The strongest and most stable vortex appears at the drag 
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maximum case; i.e., at b/e 0 = 112.5, F = 5 Hz. The cavity vortex becomes weaker 
during transition and recovers its strength during Mode III. 

A check of mass balance showed that between 93 an d 97 percent of mass 
entrance, f vdx !_„ pegion< exits the cavity. Either a slight angle of attack or 
presence of three-dimensionalities in flow are responsible for imperfections in 
the mass balance. It remains to mention that it is logical to assume that positive 
entrainment to the shear layer, at the upstream region, is the primary process 
for generating +U region and consequently causing the deflection of the shear 
layer into the cavity at the downstream corner. Relation between the positive 
entrainment and amplification of the perturbations in the initial region of the 
shear layer will be discussed in the flow model section. 

4.7. Distribution of Along y = 0 line 

k'e 

The —— distribution along y = 0 line (figures 50a through 50g) exhibits some 

interesting features. A typical distribution shows that -^—increases from zero 

at x = 0 to a peak value, and is usually followed by a stagnation flow-like 
behavior where velocity experiences a sharp decrease down to zero at x = b. 

The stagnation region often shows several local peaks and valleys. In Mode II, 
the peak value of -y— shows an increase until it reaches its maximum at 
b/e D = 112.5 , corresponding to the maximum drag case. The end of Mode II and 
the transitional region are marked by lower values for the peak of -jj— similar 

V 

to that of the - — 

Ue 


distribution. 
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Mode III does not show a substantial change in — — distribution or value of 

ijy 

peak velocity during its course (figure 51). It is important to note that -rr-lmax 
v t'a 

changes from 0.3 to 0.6 during its rise in Mode II. which is a clear indication of 

deflection of shear layer into the cavity. Consequently this deflection will cause 

the stagnation of the high velocity portion of the shear layer velocity profile at 

the downstream corner. It seems that the natural diffusion of the shear layer 

due to the cavity width increase (shear layer length), plays aless important role 

in bringing higher velocity fluid into the stagnation region. Decrease in the peak 

value of -jj - | max during the transitional region, despite the increase of cavity 
k'e 

width, supports the previous argument. In Mode 111, the averaged level of the 

~fr — I max (while having a value of 0.5) stays constant (figure 51). Diffusion of the 

shear layer into the cavity in Mode III probably plays a more important role than 
the deflection of the shear layer in bringing high velocity fluid to the stagnation 
region. 


Uv 


4.B. Evaluation of — r- Along y = 0 lane 

Uf 


uv 


The — r-term is evaluated along the y = 0 line, using the independent meas- 
Uf 


urements of ~—a.nd — — and the results are shown in figures 52a through 52g. 


U B 


u B 


A typical distribution of — s- shows a wide region of positive values for p 

Ur Ur 


uv 

r 2 
e 


which is followed by a sharp negative peak close to the downstream corner. The 

U7J 

region of positive values of — ^-contributes thrust to the net force balance while 

Ur 


USD 

the region of negative values of — =- contributes drag, (figure 53) 

Ur 
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oscillations. Even though 


-dx is always negative during Modes II and III of the 
distribution showed a fairly good mass balance 


for the cavity, the -v region was weighted by higher values of during 
integration along the y = 0 line, resulting in a positive net drag for the cavity. 


For convenience, the following conventions are introduced. The net contribu- 

|TjJ 

tion of — g-term to total drag coefficient 


= c> 


Dav 



Net contribution of 


£4? 


term to total drag coefficient 


= C, 


mw 



Total drag coefficient^ Q j. 
Therefore 


Q? - Qwr + Cjxnr 

A comparison between and 

given in the next section. 


contributions to the net drag will be 
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4.9. Shear Stress DistributionAlongy=0 line 

Figures 54a through 54g present shear stress distribution for Mode II. Shear 
stress grows along the y = 0 line, reaches a peak and vanishes steeply as the flow 
approaches the downstream stagnation point. Often two peaks appear before 


the decrease in the stagnation region starts. The peak of 


uu 
U « 


reaches the max- 


imum value of .0063 at b/do = 112.5, which corresponds to the maximum drag 
case. During the transition region the peak value of the shear stress decreases 
to 0.003. 


The start of Mode III is marked by the appearance of a wide region of high 
shear stress value close to the downstream comer. In Mode III (figure 55), the 


peak value of 


uW 
UJ 2 


increases from .007 at b/9 0 = 132 to 0.0098 at b/9 0 = 146. 


Figure 56 presents the variation of — — and — - — versus b/e 0 . It also reveals 

Q> C D 

that in Mode II, especially in the beginning stages, C ^ contributes up to 80% of 
the total drag coefficient. As b/6fe increases, the role of C a ^ becomes more 
important but its contribution never exceeds that of Mode III shows a 

dominant role for C L Immediately after transition, C D ^ constitutes 60% of 
C D and increases to 70% at the critical gap (6/0 o = 146). The message of figure 

56 is that consideration of the shear stress, -, as the only means of momen- 
ts 

turn transfer along the cavity y = 0 line will result in a 30 to 80% underestima- 
tion of the cavity drag depending on cavity length. 


Kistler and Tan (1967) made such estimations of the drag coefficient of a 
two-dimensional cavity with a turbulent upstream boundary layer. Using hot 
wire anemometry, he determined the shear stress distribution along the cavity 
mouth. His results show up to 60% underestimation of Cp with respect to Cpp. 
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It is important to recall that, due to insensitivity of hot-wire techniques to the 

direction of the flow, it is impossible to accurately estimate the term — r- along 

U e 

the y = 0 line. 

4. 10. A Possible Flow Model 

The drag coefficient of the cavity stays low as long as self-sustained oscilla- 
tions are present in the flow. Flow visualization studies have revealed that the 
shear layer (time averaged) re-attaches onto the downstream corner. Only dur- 
ing the wake mode does the shear layer stagnate below the downstream corner. 
The fact that the stagnation point (average position of maximum pressure) is 
always fixed at the downstream corner supports the flow visualization conclu- 
sion. There remains one question to answer: what mechanism keeps the flow 
from stagnating below the downstream corner and consequently results in a low 
drag value for the cavity? What follows is a possible answer to the above ques- 
tion. 


It is logical to assume that the volume of flow that is entrained by the shear 
layer is proportional to the strength of the shear forces in the shear layer. Fig- 
ure 57 reveals that the amount of fluid that is entrained by the shear layer (or 


exits the cavity \^f vdx 


+v 


follows the trend of maximum shear stress of the 


shear layer at a fixed streamwise position: e.g., x/e 0 = 65. 


On the other hand, mass balance of the cavity requires that the shear layer 
deflect into the cavity in order to supply fluid for positive entrainment to its ini- 
tial region. Figure 23 presents the location of 0.5 in the shear layer. At 

the beginning of Mode II ( b/e 0 = 83, F = 6 ), where shear stress level and growth 
is low, the shear layer displays mild curvature. At the maximum drag case 
(b/e 0 = 112.5), where the shear stress level throughout the shear layer is 
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higher, the shear layer shows steeper curvature at the initial region and a sharp 
deflection into the cavity. Therefore, the degree of the deflection of the shear 
layer depends on the nature of the instability growth in the shear layer. The 
severity of the shear layer deflection is decisive in the nature of the shear layer 
interaction with the downstream corner. 

Rockwell and Knisely (1979) show that if the respective location of the down- 
stream comer is raised by “ = 6, the rms transverse fluctuations, iT, will 

3d 

drop by a factor of 2. Sarohia (1976) demonstrated that by external mass injec- 
tion to the cavity, it is possible to deflect back the shear layer at the corner and 
weaken the strength of the cavity flow fluctuations. Based on these experiments 
it can be concluded that for given b /9 0 the displacement of the impinging shear 
layer on the corner due to deflection should occur in such a way that the feed- 
back loop would not lose its strength. Any severe external interference with the 
established self-sustained oscillations would result in an overall drop of the 
fluctuation level throughout the shear layer. These two experiments provide the 
missing link between the role of the corner and the fluctuation level. Once the 
feedback loop due to the frequency selection process is established, any inboard 
displacement of the velocity profile caused by the natural tendency of the wake 
mode will be resisted by the following mechanism. A negative displacement of 
velocity profile at the downstream corner will result in the lowering of the shear 
stress level which will be followed by a decrease in the volume of positive 
entrainment. However, a less positive entrainment would result in reduced 
deflection of the shear layer and restoration of the initial position of the shear 
layer velocity profile at the downstream corner. This mechanism will not allow 
the full velocity profile at the impinging corner to be displaced more than its 
total thickness, and as a consequence, keeps the location of C Pmnv fixed at the 
downstream corner. The strength of this restoration mechanism depends on 
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the stability of the cavity flow oscillations, which was already established by the 
frequency selection process. It is possible that because of the existence of such 
a restoration mechanism, some low frequency oscillations appear in the shear 
layer. According to the above process these frequencies need not be the har- 
monics of fundamental frequency. As stated, Knisely (1980) observed 
0.5, 0.4 and 0.6 of fundamental frequency in the shear layer over cavities. In 
the present study severe low frequency components were observed in Mode III 
(last mode) of the oscillations. Toward the end of Mode III, as a result of a 
decrease in the effectiveness of the feedback loop, the restoration force became 
weaker. Consequently, the natural wake mode of the front body overcame the 
restoration mechanism, and the shear layer had more freedom to move below 
the downstream corner and stagnate inside the cavity. 



-48- 


Chapter 5 

Control of Cavity Flow Oscillations 

One of the objectives of the present study was to develop and use a technique 
which is capable of controlling the frequency and amplitude of established flow 
oscillations. However, such a technique can be used to obtain information 
regarding the necessary conditions for the onset of self -sustained oscillations. 

This chapter will report on and discuss a series of experiments regarding the 
application of a new technique for forcing the cavity shear layer. Throughout 
these experiments, the response of the cavity shear layer, both in non- 
oscillating and oscillating modes, is studied against different forcing frequencies 
and amplitudes. 

5.1. Forcing of the Shear Layers 

Forcing of the unstable and stable modes of laminar or turbulent shear 
layers has been the topic of many previous studies. Basically, forcing techniques 
fall into two catagories: 

1) Perturbations to the bulk of either flow stream (Miksad, 1972, and Ho and 
Huang, 1982). 

2) Localized perturbations in the initial region of the shear layer (Oster and 
Wygnanski, 1982). 

In the first type of forcing, e.g,, by oscillating one of the streams, the whole 
streamwise length of the shear layer experiences forcing effects. As a result no 
accurate information can be deduced about the influence of perturbations on 
the initial regions of the shear layer. 
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Wygnanski and Oster showed that forcing effects can be localized by 
sinusoidal oscillation of a pivoted thin flap at the leading edge of the shear layer 
splitter plate. The other technique used to localize the effect of forcing is to 
introduce controlled excited (T-S)* waves, amplified by the boundary layer, into 
the shear layer. The idea of actively forcing the boundary layer was first pro- 
posed by Liepmann (1943). The most popular technique to facilitate the idea 
was the vibrating-ribbon developed by Schubauer and Skramstad (1947). The 
difficulty with either the vibrating-ribbon or the flapping edge technique for the 
present experiments is that it cannot be adapted to axisymmetric flows. 

5.2. Surface Heating Technique 

Liepmann, Brown and Nosenchuck (1982) recently showed that time- 
dependent surface heating can be used to introduce perturbations into laminar 
boundary layers. Such a technique can be easily adapted to excite T-S waves in 
an axisymmetric boundary layer such as that of the present study. The idea of 
periodic heat release into the boundary layer for the purpose of exciting T-S 
waves was also realized by Kendall (1967) in supersonic boundary layers. By 
applying periodic current to two closely spaced parallel wires flush mounted on 
the wall, Kendall was able to initiate glowing between the two wires and conse- 
quent heat release into the boundary layer. As a result of this periodic heat 
release into the boundary layer, T-S waves of up to 30 kHz were excited. 
Requirement of low pressure to attain glowing makes this technique of limited 
use for subsonic flow applications, 

For cases in water, Liepmann, et al. (1982) show that the effect of 'heating 
and cooling of the surface is roughly equivalent to suction and blowing or nega- 
tive and positive surface displacements," with 


C 


ToUnnen-SchHchting 
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where fi is dynamic viscosity, T is the temperature of the fluid, Pr is the 
Prandtl number, w is the width of the strip, and R x is the Reynolds number. 
Such an equivalent normal surface oscillation would periodically redistribute 
vorticity. Depending on the state of the boundary layer, the external perturba- 
tions would either be dampened or amplified exponentially. Linear theory can 

predict the exponential growth of the (T-S) waves, based on the parameters -jj— 
and Red* for a given pressure distribution. 

5.2.1. Application of T-S Waves as a Means of Forcing of the Shear Layer. 
Forcing T-S waves in boundary layers has been practiced for a long time. How- 
ever, to our knowledge, no use has been made of T-S waves to force a separated 
shear layer resulting from such a boundary layer. In the present study, 
inspired by the experiment of Liepmann, et al., periodic heating of a thin film 
was used to excite the T-S waves in the boundary layer upstream of the cavity 
shear layer, benefitting from amplification of the waves before their introduc- 
tion into the shear layer (figure 58). One can imagine that the effect of the T-S 
waves or v 1 component at the upstream corner of the cavity would be similar to 
that of the oscillating flap at the leading edge of the splitter plate in Oster and 
Wygnanski’s experiment. It will be seen later that the frequencies needed to 
excite the free shear layer are available in the range of boundary layer T-S 
waves. Thus, a transformation of boundary layer T-S waves to the shear layer 
Kelvin-Helmholtz waves is possible. 
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5 . 2 . 2 . Strip Heater Set-up. Initially three strip heaters were mounted on the 
ellipsoidal nose of the model (figure 59) at s = 9.2 cm, 12 cm, and 15.5 cm where 
s is the distance on the model surface from the nose of the model. Double stick 
scotch tape was used to electrically isolate heaters from the model. The first 
and third heaters were stainless steel ribbons 3.7 mm (.125 inch) wide, .051 mm 
(.002 inch) thick, wrapped around the front nose of the model. The second 
heater was narrower, .08 cm (.2 inch) wide. Each heater entered the model from 
a flush mounted round adaptor, designed to connect the heaters to the current 
supply lines. The resistance of the strip heaters was 2-4 ohms. 

A high quality commercial dual channel power amplifier (Hafler DH-500) was 
used to supply current to the strip heaters. DH-500 is a 500 watts per channel 
amplifier which can be easily modified to provide frequencies even below 1 Hz 
while delivering adequate power to the strip heater. An HP Function Generator 
(3325A) was used to provide the selected frequency to the power amplifier. It 
should be mentioned that due to the quadratic dependence of Joule heating on 
the voltage, the input frequency to the power amplifier had to be half of the 
desired forcing frequency. 

5 . 2 . 3 . Important Parameters of Strip Heater Technique. Several parameters 
are important for the efficiency of a strip heater in forcing both the T-S waves 

W 

and the shear layer. Laepmann, et al. (1982) stated - 7 - (where W is the strip 

0 

heater width and 6* is the boundary layer displacement thickness) is an impor- 
tant parameter in effectiveness of strip heater operation on the local forcing of 
T-S waves. The state of the free stream pressure gradient is another important 

W 

factor. In the present study, the width-to-boundary layer thickness ratio, — 

0 

was 5.6 for the first strip heater and maintained at 3 for the second and third 
heaters. 
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The preliminary experiments indicated that for a fixed input frequency and 
power level, the first heater, s = 9.2 cm, was far more effective in forcing the 
free shear layer than the second and third heaters, set at s =12 cm and s = 15.5 

W 

cm, respectively. Also, it was noticed that having the same ratio the 

<5 

second heater operated more effectively than the third heater. Figure 59 reveals 
that (T-S) waves initiated by the first and second heaters benefitted from a posi- 
tive pressure gradient as well as a longer period of amplification than that of the 
third heater. 

In the final version of the experimental model, a single strip heater was flush 

W 

mounted at s = 9.2 with 5.3, 

6 

5.3. External Forcing of the Non-Oscillating Cavity Flow 

The receptivity of the cavity shear layer to externally imposed disturbances 
was studied. Two different cavity widths were examined, with corresponding 
b/e = 66< 6/0 1 min and b/e- 77 < 6/01,^. Throughout this section these 
selected widths will be called case A and B respectively. In both cases edge velo- 
city was fixed at U B = 22cm. /Sec . Cavity flow response was measured in terms of 
the maximum shear stress that was developed in the shear layer at a station 
,254 cm (.1 inch) upstream of the downstream corner. 

Spectral analysis of the response velocity fluctuations indicated a one to one 
relation between forcing and the response frequency. For Case A, figure 60 
presents the response of the cavity shear layer to a wide band of the forcing fre- 
quencies at different levels of forcing power. The results indicate that the level 
of shear layer response increases as the forcing power increases. For forcing 
levels above 9 watts, resonance peaks appear at the forcing frequency of 7.2 Hz. 
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An overall phase measurement revealed that Case A satisfies the phase cri- 
terion 2, at the resonance (F = 7.2 Hz), thus verifying it as the 

2tt U c 

second mode of oscillation. The absence of forcing caused the resonance mode 
to disappear, thus revealing that the flow oscillations are not self-sustaining. 
Figure 60 also reveals that the existence of resonance depends on a minimum 
amplitude threshold. 

Figure 61 presents the response of the cavity shear layer for Case B at a fixed 
forcing level of 9 watts. This case shows a resonance peak at the forcing fre- 
quency of 6.2 Hz with g~ = 2- Appearance of resonance peak for Case B 

at a lower forcing level than that of case A indicates that the threshold level 
decreases as b/e 0 increases. In Case B, the effect of initial lower amplitude of 
the forcing was compensated by longer cavity length in such a way that reduced 
the required forcing level for the resonance condition. The effect of higher forc- 
ing power on the exponential growth of the perturbations was to increase the 
general level of fluctuation growth (figure 62a). Figure 62b presents the growth 
behavior of the different frequencies for a fixed forcing level (9 watts) for Case 
B. 


It can be concluded that external frequencies can be imposed on naturally 
non-oscillating cavity flows, within the receptivity range of the shear layer. 
These frequencies can cause the cavity flow to resonate if they satisfy the follow- 
ing conditions: 

1) 7 ^— = N for the given width (b); and 
2rr 

2) A > Afh where A is the amplitude of external frequency at the down- 
stream corner and Afh, is the threshold amplitude. It is logical to propose that 
as the cavity width increases, the threshold level decreases to such an extent 
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that a flow background frequency, which satisfies the phase criterion and has 
sufficient amplitude, will initiate the self -sustained oscillation. 


It is now clear that due to the absence of proper frequencies with sufficient 
intensities in the flow background, the first mode did not appear naturally in 
this study. However, by selecting a proper forcing frequency and initial ampli- 


tude, the first mode of oscillations 


-Z-=££-» i 

2rr U c 


was simulated for Case A. 


This was accomplished at a forcing frequency of 4 Hz, a power level of 35 watts. 
A flow visualization of this case is shown in figure 63b. In comparison to figure 
63a, which is the same case but unforced, two features are noteworthy. The first 
is the appearance of a small amplitude wave in the shear layer of the forced 
case. The second is the interaction of this wave with the downstream corner and 
the resultant strong recirculating flow inside the cavity. 


It is interesting to note that the response of the cavity flow to external forc- 
ing is not instantaneous. Figures 64a and 64b present simultaneous traces of 

X X 

the velocity fluctuations taken at stations corresponding to — = 10 and — -= 60, 

0 0 

respectively. The cavity was in the non-oscillating mode with 6/9 = 66. Figure 
62a indicates that for an external forcing frequency of 7 Hz, simulating Mode II, 
it takes 3 cycles of the oscillation before the peak amplitude establishes itself. 
When the forcing is turned off it takes many more cycles for the amplitude of 
oscillations to diminish. 


5 . 3 . 1 . Effect of External Forcing on Cavity Drag in the Non- Oscillating 
Mode. Case B, at the resonance, had a drag coefficient of 0.0012. This drag 
coefficient was one order of magnitude higher than that of the same case but in 


the absence of forcing, Distributions of 


uv 

U* 


and 


uv 

U* 


along the y = 0 line for 


both the forced and unforced cavity are given in figures 65a and 65b. These 
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figures reveal that the difference in the drag coefficient of forced and unforced 
cases can be contributed to the development of a strong shear stress level in the 
forced case. 


5.4. Effect of External Forcing' on a Naturally Osc illatin g Cavity Flow \b/e > b/e]^ 

Once the cavity flow established the self-sustained oscillations 
> b /&q | mjjj , the flow resisted any external disturbances. Thus in order to 

interact with the flow oscillations, it was expected that the amplitude level of 
external forcing must be comparable to that of the self-sustained fluctuations 
in the flow. 


A non-dimensional cavity width of b /9 0 = 95 was selected for an external 
forcing experiment. The cavity flow was oscillating at 5.4 Hz in Mode II. An 
external forcing frequency of 6 Hz was introduced to the shear layer. Depicted 
in figures 66a through 66g are the spectral distributions of velocity fluctuations 
for several levels of forcing as the forcing power was changed from 2 to 27 watts. 

Spectra reveal that as the amplitude of the forcing increases, the amplitude 
of natural oscillation decreases and eventually disappears. At the forcing power 
of 40 watts, the forced frequency becomes the dominant frequency of the oscil- 
lations. Oscilloscope traces of the signal for the case when both natural and 
forced frequency have the same amplitude, as seen in figure 66d, are given in 
figure 66h. Strong modulation of the signal indicates that both frequencies are 
present simultaneously. Note that in the co-existence region, figures 66b 
through 68f, the peak amplitude of both frequencies is lower than that of the 
single frequency regions of figures 66f and 66g, 

Miksad (1973) in his study of dual excitation of a free shear layer observed a 
similar interaction of two frequencies. He reported that if the excitation 
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amplitudes were adjusted so that both disturbances reached finite amplitude 
together, the fined amplitudes of both components tended to be reduced, but 
not as strongly as when, one component was clearly dominant. Miksad concluded 
that the distortion of the mean flow by the stronger frequency made velocity 
field unsuitable for the growth of the weaker frequency. Similar behavior had 
also been observed by Sato (1970) in symmetric wakes. 

Important to our study is the fact that the amplitude of the forced frequency 
in figure 66g, where natural frequency just disappears, is equal to that of the 
natural frequency in the absence of the forced frequency (figure 66a). This fact 
establishes the amplitude of the natural frequency as the level which should be 
exceeded by the amplitude of forced frequency to change the frequency of the 
oscillations. With this rule it is possible to impose forced frequencies on the cav- 
ity shear layer within the receptivity range of the shear layer. These forced 
oscillations are not, in general, self-sustaining. During the mode switching 
region (hysteresis region), self-sustaining oscillation in a desired mode could be 
achieved by providing a proper frequency within the hysteresis region of that 
mode. Thus, the mode switching phenomena could be eliminated by external 
forcing. Transition to the wake mode could be prevented by enhancing the 
amplitude of natural frequency through external forcing. It is possible to 
extend Mode III, last mode in this study, by a 20 momentum thickness through 
external forcing. However, once the cavity entered the wake mode, it was not 
possible to reestablish cavity flow oscillations by external forcing. 

A map of possible interaction with the cavity flow oscillations observed in this 
experiment is given in figure 67. 
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5.5. Cancellation Experiment 

Cavity width, of b/e^ = 82, with frequency of 6 Hz, was selected for the cancel- 
lation experiment. The signal from a downstream location was fed to a function 
generator (Tektronix FG504) which had a built-in phase lock loop. It was possi- 
ble to lock to lower harmonics of the input signal. Phase variation of the output 
signal with respect to the first subharmonic of the input signal was possible 
through an independent manual adjustment knob. The power amplifier was then 
supplied by the output of the function generator. The experimental set up is 
sketched in figure 68. 

The following procedure was conducted in the cancellation experiment. 

1) Output amplitude of function generator was set at 0. 

2) Phase lock loop was locked to first subharmonic of the input signal . 

3) Knowing the threshold amplitude, the power was fixed at a level 
which was enough to simulate conditions as depicted in figure 66d. 


4) Step (3) was repeated for different phase angles between input and 
output of the function generator (phase lock loop) . 


Due to the small rapid phase changes of the cavity flow oscillations it was not 
possible to have a phase lock situation for an extended period of time. Occa- 
sionally a correct phase angle could be found where the forced and natural fre- 
quencies would interact in such a way that wide band rms fluctuations were 
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reduced. A sample of such a cancellation is shown in figure 69. A reduction in 
rms fluctuation was achieved by a factor of 2 during 30 cycles of the cavity 
shear layer oscillations, Liepmann and Nosenchuck (1982) report a factor of 2-3 
in the reduction of rms fluctuations of the natural (T-S) waves in the flat plate 
boundary layer. 
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Chapter 6 

CONCLUSIONS 


6.1. Summary of Results 

(a) Investigation of the flow field in the oscillating mode reveals that the 
momentum thickness of the cavity shear layer grows linearly with x/e 0 . The 
linear growth of the momentum thickness may indicate the existence of self- 
similarity in the flow. However, the self-similarity argument for the free shear 
layer requires a constant maximum shear stress along the layer. In contrast, 
cavity shear layers show an exponential growth of the maximum shear stress 
with x/Qq. The cavity shear layer in the non-oscillating mode does not develop 
any noticeable shear stress along its span. The velocity measurements also ver- 
ify the existence of a recirculating flow inside the cavity. Streamwise and 
transverse velocity fluctuations as high as those of the shear layer are detected 
in the recirculating region. 

(b) Upstream laminar boundary layer has a shape factor value of 2.5. The 
effect of the cavity on this profile is found to decrease the shape factor to a 
value of 1.54, which is close to that of the turbulent boundary layer. 

(c) The drag measurement results presented in figure 39 indicate that the 
drag coefficient of an axisymmetric cavity with laminar upstream boundary 
layer remains low as long as the cavity flow stays in a non-oscillating or self- 
sustained oscillating mode. Results also reveal that the drag of the cavity in the 
non-oscillating region is less than the case in which the cavity was replaced by 
a solid surface. 
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(d) Flow visualization study (figure 16) in conjunction with the drag measure- 
ments reveals that for a long enough cavity width, b/R> 1.2, cavity flow exhibits 
a wake-type behavior. Also, the average drag coefficient in this mode, Co » 0.35, 
is one order of magnitude higher than the highest drag coefficient in the oscil- 
lating mode, Co = .036, 

(e) The most important result is presented in figure 40. The drag coefficient 
of the cavity exhibits an exponential dependence with b/e^ in Mode II of the 
oscillations followed by a decrease in the same mode and another course of 
exponential growth in Mode III. An independent estimate of the drag coefficient 
through measurement of the shear stress and mean momentum transfer terms 
Uv along the y = 0 line confirmed the exponential growth of the cavity drag. 
Furthermore, these measurements revealed that the uv term is responsible for 
30% to 80% of the net drag coefficient of the cavity. 

(f) Study of the cavity pressure distribution revealed that the pressure rise 
close to the outer region of the downstream comer is responsible for the drag 
coefficient of the cavity. 

(g) In the oscillating modes of the cavity flow, the maximum pressure 
coefficient ( C Pmia ) always occurs at the upstream corner and stays so until the 
wake mode appears. In the wake mode the position of maximum pressure 
coefficient (^ max ) moves toward the center of the cavity as b/6o increases. 

(h) The existence of the feedback loop introduces an extra stability to the 

cavity flow field which includes a fixation feature of at the upstream 

comer and a resultant low drag for the cavity. 

(i) The distribution of the shear stress along a straight line connecting two 
comers of the cavity is shown to have direct effect on the positive transverse 
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velocity (-H7) distribution on the same line. It is also shown that the deflection 
of the shear into the cavity is related to the distribution of positive transverse 
velocity through the mass balance of the cavity. 

(j) Stronger deflection of the shear layer into the cavity causes a higher wu 
region to develop close to the downstream corner and the drag of the cavity to 
increase. Thus, shear stress development along the shear layer has a direct 
effect on the drag of the cavity as well as an indirect effect through its role on 
the development of the wv term. 

The type of shear layer interaction with the downstream comer which 
depends on the shear layer deflection is decisive in the development of the shear 
stress throughout the cavity shear layer. Based on the findings of previous 
experiments, it is proposed that inward displacements of the shear layer with 
respect to the downstream corner will result in an overall drop of the shear 
stress level. As a result of this overall shear stress drop, a less positive region of 
transverse velocity develops. Consequently, the mass balance of the cavity 
requires less deflection of the shear into the cavity. It also requires that the 
shear layer deflect back and restore the position of the shear layer. This 
mechanism keeps the position of Cp „„ fixed at the upstream corner and it 
might also introduce some low frequency components to the oscillations. 

(k) The control of the cavity flow oscillation is facilitated by external forcing 
of the cavity shear layer. It is shown that the excited T-S waves can be used to 
force the shear layer generated from separation of respective boundary layers. 
Thus, a transformation of the T-S waves of the boundary layer to the Kelvin- 
Helmholtz waves of the shear layer is possible. 

(l) For short cavities, i.e., non-oscillating mode b/9o 82, periodic flow oscil- 
lation at a wide range of frequencies was obtained. A resonance condition was 
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ob served for a selected forcing frequency if it 

1) satifled the phase criterion -N 

2tt 

and 

2) its amplitude exceeded a certain threshold amplitude. 

For a fixed amplitude the threshold level increases as b decreases. The forced 
resonance cases are not self-sustaining. 

(m) Previous conclusions indicate that, for a given b/e 0 , the first frequency 
in the flow background that satisfies the phase criterion and has sufficiently 
high amplitude will start the self-sustained oscillations. Thus the determination 
of the b/e olmin strongly depends on the turbulence content of the flow back- 
ground. 

(n) Cavity drag in the non-oscillating mode can be increased by one order of 
magnitude through external forcing. Non-oscillating cavity flow at b/e a = 77 
produces a drag coefficient value of 1.3 x 10 -4 . External forcing at 6 Hz and 
power level of 9 watts increases the drag coefficient to a value of 1.2 x 1CT 3 . It is 
also possible to generate shear stress levels several orders of magnitude higher 
than that of the laminar shear layer by external forcing. 

(o) For naturally occurring oscillations it is possible for two waves to co-exist 
in the shear layer (natural and forced). Spectral analysis, however, shows that 
the amplitude of the forcing frequency has a direct effect on the growth of the 
natural frequency. External forcing made possible the ability to completely sub- 
sidize the natural self-sustained frequency by a non-self-sustained frequency. It 
is shown that the amplitude of the established self-sustained oscillation is the 
amplitude level that should be exceeded by the amplitude of the forced fre- 
quency for complete switching of self-sustained to forced frequency. 
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(p) Mode switc hing or transition to the wake mode can be delayed by enhanc- 
ing the self-sustained oscillations through external forcing. Once the wake 
mode appears, the self-sustained oscillations cannot be reestablished. Different 
mode extension and simulation is possible by providing the proper forcing fre- 
quency to the cavity shear layer. 

(q) The most important result of external forcing is the possibility of cancel- 
ling or damping the established self-sustained oscillations. An external forcing 
frequency equivalent to that of the cavity with an amplitude less than threshold 
level which has a proper phase angle respective to that of the self-sustained one 
was able to reduce rms fluctuations by a factor of 2. 

6.2. Concluding Remarks 

This study has disclosed that the presence of self-sustained oscillations helps 
in keeping a low level of drag coefficient for the cavity. A more detailed study of 
the role of the downstream corner on the overall development of the shear 
stress along the shear layer and on the drag of the cavity is desirable. It is 
important to clarify the effect of the cavity presence on the net drag coefficient 
of the system. 

With respect to the active control of the cavity flow oscillations, if a more ela- 
borate feed-back system is used, a high percentage of fluctuation cancellation 
might be attainable. 



- 64 - 


Appendix A 

Determination of Overall Phase Difference for the Cavity 

For the purpose of obtaining phase information between two time dependent 
signals, the cross-correlation function between these two signals should be car- 
ried out. By cross-correlating two signals, the relative non-dimensional phase 
change between the location of the reference probe and the second probe can be 
determined. This is done, as shown in sketch (A-l), by measuring the time 
shift, along the positive T-axis, between two corresponding peaks in the cross- 
correlation function of two signals and dividing this time by the average period 
between peaks. 



Then the non-dimensional phase difference between two locations separated by 
Ax will be 
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A <p 
2n 


r SHIFT 

r AVERAGE PERIOD 


( 1 ) 


The wavelength of the oscillation can be determined as 


\ = Az • 


Srr 
A <p 


( 2 ) 


Knowing the fundamental frequency of oscillation, the phase velocity can be cal- 
culated from 


U 0 = f\ 


( 3 ) 


Overall phase difference between the two corners is given by 


9 = 


b 

Ax * 


A 93 


( 4 ) 
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Appendix B 

Laser Doppler Yelocimeter Optical Arrangement 

A two watt argon-ion laser is used as the light source. Past a collimating lens, 
a color separator provides two lines of the laser, the 514 nm green line and 488 
nm blue line. Each beam is split into two pairs of equal intensity beams by a 
polarization rotator and a beam splitter. Each pair passes through its desig- 
nated frequency shifter. Non-desired diffracted beams are blocked by two slit 
mechanisms. The measuring volume is formed by focussing the four parallel 
beams by a single 598 mm focal length lens. The collimator is used to minimize 
the lateral extent of the focal volume. 

The scattered light from the particles is collected by a 600 mm focal length 
lens. A dichroic mirror is used for color separation, with a focussing lens used 
to focus the scattered light of each of the two colors into a pin-hole (250 (j, m) 
and photomultiplier tube assembly. A first-surface flat mirror is used to rotate 
the beam 90°, A second mirror is used to deflect the scattered light into the col- 
lecting optics. 

The component of the velocity vector measured by the LDV and the Doppler 
constant, the conversion factor from frequency to velocity, is determined by the 
transmitting optics. The green beams' pair enters the test section in a plane 
forming an angle of 0° with the downstream direction, while the blue beams' 
plane forms an angle of 90°. Therefore, streamwise velocity component is meas- 
ured by the fringe system formed by the green beams while velocity components 
perpendicular to the mean free stream are measured by the fringe system 
formed by the blue beams. Figure B-l presents the schematic of focal volume 
formation. The Doppler constant which is related to the fringe spacing within 




/ 
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the focal volume is given by the equation 


ic = 


X 


2 sin 


a 

2 


where X is the light wavelength and a the angle at which the beams cross in the 
focal volume. The angle a = 4.72° was used which gives Doppler constants for 
the green and blue components kq = 6.24 (m/Sec) ( MHz)~ l and 

kb - 5.92 (m/&ec) respectively. 
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Appendix C 

Velocity Bias and Bra gg Cell 

Fringes of the measuring volume are formed by the interfering effect of the 
two beams crossing. These fringes are stationary and do not produce Doppler 
frequency for a stationary particle in the measuring volume. 

Shifting the frequency of one beam can produce a system of moving fringes. 
Then, a stationary particle in the measuring volume will give an output fre- 
quency from the photodetector. If the velocity increases in the direction of the 
fringe movement, the frequency will decrease. If it is opposite to the fringe 
movement, the frequency will increase. Therefore knowing the direction of the 
fringe system movement will discern the particle direction. 

In this experiment, two acousto-optic cells (Bragg cell) were used to fre- 
quency shift one beam of each color by a fixed 40 MHz. Each cell is installed in a 
cylindrical housing that aligns with the optical axis. A compensating wedge is 
used to correct for shifted beam deflection, and a tilt adjustment allows shifting 
light intensity to the deflected beam, as much as 80%. The cell is derived by a 40 
MHz, 2 watt signal provided by a power amplifier. A frequency of 40 MHz is often 
too large. A downmix circuit permits decrease of the effective frequency shift of 
the photodetector output to values from 2kHz to 10 MHz in a 2, 5, 10, sequence. 
Bragg shifts of nominally 30 and 50 kHz were used. In this study the green (U 
component) channel was shifted against the mean flow, while the blue channel 
(V component) was shifted in the positive y direction (figure B-l). 



-69- 


Appendix D 

LDV S ign al Processing 

The burst processing encompasses the required electronics for the actual 
measurement of the particle velocity. A particle burst is considered valid if it 
satisfies the following conditions: 

1) Its amplitude crosses zero and exceeds the threshold level (about 50 mV) . 

2) For zero crossing (with hysteresis), the signal should be more than 10 mV. 


The pulse that occurs on a zero crossing after the threshold has been exceeded 
is the signal of the interest. At this point the leading edge of a time envelope is 
marked (by the processor) and gates a high speed clock (250 MHz) into a 
counter. The length of the envelope is determined by the total number of cycles 
exceeding the threshold. The first cycle that exceeds the hysteresis level but 
not the threshold level determines the trailing edge of the envelope. The trailing 
edge stops the counter and the clock. The input frequency is then: 


INPUT FREQUENCY («) = 


ENVELOPE CYCLE NUMBERS 
TOTAL ENVELOPE TIME 


The resolution of time measurement is 1 nsec. In this experiment only bursts 
with more than 16 cycles were used for processing. A schematic of valid burst 
processing is given in figure D-L The output of the burst processor is in digital 
form. It feeds into the computer interface that incorporates the coincidence 
test between the two processors. The coincidence test consists of the initiation 
of a time window after completion of processing by one of the processors. If 
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data from, the other processor are not received before the end of the window, the 
data are discharged; otherwise coincidence is attained and the data from both 
processors are transferred to the computer. 

A minicomputer was used to acquire and store the data. A computer program 
was developed for data acquisition using the LDY system. This computer pro- 
gram accomplishes two tasks. It generates a header which contains information 
relevant to the conditions of the experiment and the LDY setup. It also acquires 
a buffer of 16,000 words of data. This buffer size corresponds to a total of 3200 
velocity measurements in two-component measurements, since two words are 
generated for each component and one word was used for the time interval 
information between data points. 
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Appendix E 
LDV Data Processing 

The following expressions cire used to calculate, the mean, second moment 
and cross-correLation between two veLocity components. 
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where At* is the time between data points and a is the frequency. The sub- 
scripts G and B refer to the frequencies measured by the green and blue 
beams, respectively. The equations relating the mean and the correlations of the 
various velocity components are 


u — 


~ k g , v - u b ~ k b 


(E— 6) 


u 2 = (w 2 c - 2 u FS C + F5 2 c ) k? g (E-7) 

g — 2 S3 5 FSg + /’S'jj K Z g (E-8) 

WU = jwcCJp — CJ G FSg — UgFSc + Kg X Kq (E-9) 


where FSc and FSg are frequency shifts of green and blue channels, 
respectively, and k. is the frequency to velocity conversion factor. 

Rms fluctuations and shear stress values can be obtained using the following 
relationships: 

iT = - u z (E-10) 

= y/v 2 — v* 


v 


(E-ll) 



No corrections were made to the data presented in this study. 
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Fig. 1 Water tunnel schematic 
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Fig. 2 Power spectrum of the velocity fluctuations 

(a) In the absence of the model. 

(b) Typical cavity flow oscillation. 
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Fip. 3 Schematic of the experimental model. 
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Fig. 4 Schematic of the pressure sensing set-up and operational 
table of the valves (Table 1). 
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Fig. 5 Photograph of the laser Doppler velocimeter set-up and the test section. 
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Transmitting Optics Number of fringes, Nf r = 48 

Fig. 6 Relative location of the LDV with respect to the model and characteristics of the 
measuring volume (Table 2). 
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Fig. 7 Frequency variation of cavity shear layer with varying cavity width 
constant free stream velocity. 
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8 Velocity fluctuation spectra for increasing 
cavity width from Mode II to Mode III. 
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(a) Mode II, F = 6 Hz, vertical scale .02 volt/div 



(b) Transition F = 4.6 Hz Mode II 

F = 6.4 Hz Mode III, Vertical scale .02 volt/div 



(c) Mode III F = 6 Hz vertical scale .04 volt/div 


Fig. 9 Oscilloscope traces of cavity shear layer oscillations. 
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Fig. 12 Overall phase difference between two corners of the 
cavity for Modes II and III of the cavity flow 
oscillations. 
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Fig. 13 Flow visualization of the cavity flow in non-oscillating 
mode; b/0o = 66. 
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(b) Modem b/X = Yp- = 3 

Uc 

Fig. 14 Flow visualization of the cavity shear layer; flow is 
frozen at 1/1000 sec. 
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Fig . 15 Frequency variation of the cavity shear layer with — , for different 
cavity geometries and free stream velocities. 
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Fig. 18 Non-dimensional mean streamwise velocity profiles; b/0 n = 103, 
F = 5.1 Hz, Mode II. 
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Fig. 19 Non-dimensional mean streamwise velocity profiles; b/0 o = 130, 
F = 6 Hz, Mode III. 


Vo =30 Vfl =65 

°o © | v o 


96 



Fig. 20 Non-dimensional mean transverse velocity profiles; b/0 
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Fig. 21 Non-dimensional mean transverse velocity profiles; b/0 o = 103, 
Mode II. 
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Fig. 24 Locus of points (y 0 . 5 ) where u/U„ = 0.5 for cavity widths 0 b/0 o = 85 and 
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Fig. 26 Shear stress profiles; b/9 
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Fig. 27 Shear stress profiles; b/9 0 103, F = 5.1 Hz, Mode II. 
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Fig. 29 Growth of the maximum shear stress in the cavity shear layer 
with downstream distance; 

Mode II , |p 
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Fig. 30 Growth of the maximum shear stress in the cavity shear layer 
with downstream distance; 

Fb 

Mode III, -jp - 3. 
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Fig. 32 Growth of the maximum streamwise and transverse rms velocity fluctuations in the cavity 
shear layer with downstream distance; 

Mode III, — 3 . 
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Fig. 34 Non-dimensional transverse velocity profiles of the cavity flow; b/fi = 103, 
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Fig. 36 Non-dimensional transverse velocity fluctuation profiles of the cavity flow; 

Fb 

b/0 o = 103, F = 5.1 Hz, Mode II - 2. 
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Fig. 37 Shear stress profiles of the cavity flow; b/0 = 103, 
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Fig. 38 Effect of the downstream corner on the velocity fluctuations 
and shear stress profiles of impinging cavity shear layer; 

b/6 0 = 103, F = 5.1 Hz, Mode II - 2. 
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Fig* 39 Variation of the ax i symmetric cavity drag coefficient with non-dimensional 
cavity width, b/0 o . 
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Fig. 40 Exponential dependence of the cavity drag coefficient on b/0 o . ® '-'pp * A . 

♦ ~ estimated C based on boundary layer friction in the absence of the cavity. 



Fig. 41 Pressure coefficient distribution in Mode II of the 
oscillations 

• Upstream face pressure coefficient distribution. 
©Downstream face pressure coefficient distribution. 
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Fig. 42 Pressure coefficient distribution in Mode III of the 
oscillations 

• Upstream face pressure coefficient distribution. 
O Downstream face pressure coefficient distribution 
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Fig. 46 Variation of pressure coefficient at the middle of the downstream face with b/9 



Fig. 47 Variation of the cavity drag coefficient with b/R; © present study; • Koenig (1978) 
optimum drag coefficient of disk-cylinder combination* 
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Fig. 48 Continued. 
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Fig. 52 Distribution of the 
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Fig. 52 Continued. 
















contribution to the net cavity drag 




Fig. 57 (a) 
(b) 
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Fig. 59 Pressure distribution over the ellipsoidal nose cone and 
respective location of the strip heaters. 
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Fig. 61 Variation of the maximum shear stress in the shear layer 
with forcing frequency, constant power. 
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Fig. 63 Flow visualization of a naturally non-oscillating cavity flow, b/9 - 66: 

(a) unforced, (b) forced at F ■ 4 Hz, simulating Mode I of the oscillations. 
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Fig. 64 Oscilloscope traces of forced fluctuations: 
(a) turn on; (b) turn off; 

b/6 0 = 66, F = 7 Hz, 

c 


^ 2 . 




along the y = 0 line i 
♦ forced at 6 Hz. 
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Fig. 66 Velocity fluctuation spectra, presenting interaction of the 
natural and forced frequencies in the cavity shear layer for 
different power levels (P^ ) » (h) scope trace of Case (d) . 
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Uc 
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Fig. 67 Possible interactions with the cavity flow through 
external forcing 

1) Mapping of Mode II into Regions 1 and 2. 

2) Mapping of Mode III into Region 3. 

3) Extension of Modes II and III. 

4) Delay of the wake Mode. 

5) Simulation of Mode I. 
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Fig. 68 Schematic of the cancellation set-up. 



Fig. 69 Damping of natural oscillations through 
external forcing 

(a) natural oscillation 

(b) strip-heater on. 
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Fig. B-l Measuring volume formation. 
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Fig. D-l Valid burst processing and coincidence 
window operation. 




